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For  untended  manufacturing  operations,  in-process 
monitoring  of  tool  breakage  plays  a critically  important 
role.  A system  for  sensing  milling  cutter  breakage  has  been 
developed  using  a deterministic  approach.  It  is  based  on 
recognizing  the  changes  of  cutting  force  patterns.  The 
regular  periodic  variation  per  tooth  is  filtered  out  by 
synchronized  sampling,  then  averaging  per  tooth.  A 
distinct  tool  breakage  feature  is  clearly  available  after 
differencing  the  average  cutting  force  per  tooth. 

Due  to  practical  difficulties  for  cutting  force 
measurements  in  milling,  the  use  of  vibration  signals  is 
explored.  It  has  been  found  that  the  displacement  signal 
of  the  spindle  is  suitable  to  replace  the  cutting  force 
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signal  for  the  purpose  of  cutter  breakage  detection  using 
this  algorithm. 

Numerous  cutting  tests  have  been  done.  It  has  been 
shown  that  the  displacement  signal  of  the  spindle  is 
effective  in  the  recognition  of  cutter  breakage  under 
varying  cutting  conditions.  Furthermore,  this  processing 
method  is  also  useful  for  high  speed  milling  and  can  be 
easily  implemented  in  real-time.  Finally,  the  threshold  for 
cutter  breakage  detection  is  further  investigated  to  make 
it  independent  of  changes  in  cutting  parameters. 
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CHAPTER  1 
INTRODUCTION 


Unmanned  machining  centers  (UMC)  have  become  a vision 
for  the  future  of  manufacturing  technology.  This  could 
tremendously  save  manpower  and  maximize  machine  uptimes  so 
as  to  increase  productivity.  However,  the  feasibility  of 
using  the  unmanned  machining  centers  is  greatly  dependent 
on  the  development  of  diagnostic  systems  for  machining 
centers.  One  of  the  most  important  diagnostic  systems  is 
the  sensor  for  detection  of  tool  failure.  Generally 
speaking,  there  are  three  types  of  tool  diagnosis  [1]. 

The  first  type  is  the  tool  wear  sensor.  This  sensor 
should  be  sensitive  to  the  amount  of  tool  wear.  Therefore, 
it  would  be  possible  to  calculate  the  rate  of  tool  wear. 
The  cutting  tool  could  be  removed  from  use  before  it 
reaches  enough  wear  to  produce  excessive  cutting  force  or 
temperature,  which  would  eventually  lead  to  complete 
failure  of  the  cutting  edge. 

The  second  type  is  the  sensor  for  immediate  detection 
of  any  tool  failure,  such  as  chipping,  cracking  or  massive 
failure  of  a cutting  edge.  This  sensor  would  provide  an 
emergency  action  quickly  stopping  the  machining  operation 


1 


2 


in  order  to  avoid  further  damage  to  the  tool,  workpiece,  or 
tool  holder. 

The  third  type  is  the  protection  device.  This  device 
has  no  direct  information  about  the  cutting  tool.  It  is 
designed  to  protect  the  machine  or  the  cutting  tool  from 
damage  due  to  the  overloads  regardless  of  the  cause. 

Milling  is  an  extremely  important  machining 
operation.  However,  an  efficient  and  reliable  tool 
diagnostic  system  for  milling  is  still  not  available.  In 
milling  operations,  the  development  of  tool  wear  is 
usually  gradual,  not  sudden.  Hence,  the  monitoring  of  tool 
wear  could  be  considered  using  off-line  techniques  [2]. 
Vision  and  pattern  recognition  systems  [3]  might  be  used  to 
repeatedly  observe  the  edges  of  the  milling  cutter  between 
cutting  operations  in  order  to  assess  the  tool  wear  width. 
If  the  tool  wear  reaches  a certain  level,  a precautionary 
action  could  be  taken  so  as  to  avoid  any  further  damage. 

However,  the  occurrence  of  tool  failure  is  unavoidable 
for  milling  with  heavier  chip  loads  or  with  workpiece 
materials,  such  as  steel,  that  are  hard  on  the  cutting 
edge.  If  tool  breakage  or  chipping  occurs  suddenly  by  any 
cause,  the  damaged  cutting  edge  is  still  in  use.  The  local 
breakage  is  accumulated  and  then  spread  to  the  following 
teeth  as  a chain  reaction  because  of  the  overload  of  the 
following  teeth.  As  a result,  catastrophic  consequences  for 
the  tool,  the  workpiece,  and  the  machine  might  happen. 
Therefore,  tool  failure  sensing  is  much  more  important  than 


3 


tool  wear  sensing,  especially  for  high  speed  high  power 
milling.  The  task  of  this  research  is  therefore  to  attempt 
to  achieve  the  second  type  of  tool  diagnosis  as  mentioned 
above,  i.e.  immediate  detection  of  tool  failure  in  milling. 

Under  these  considerations,  the  requirements  of  a tool 
breakage  sensor  system  [4]  can  be  summarized  as  following: 

1.  Reliability,  on  the  safe  side: 

100%  alarm  on  tool  damage  beyond  a set  level. 

2.  Robustness: 

High  sensitivity  to  breakage,  very  low  sensitivity 
to  transients  (entry,  exit,  interrupted  cuttings) , 
and  to  variations  in  cutting  parameters  (workpiece 
material,  feed,  speed,  radial  immersion,  up/down 
milling) . 

3.  Highly  automated  threshold  setting: 

Adjustable  to  formulate  and  set  an  alarm  threshold 
under  varying  cutting  conditions. 

4.  Practicality  of  sensor: 

Easy  to  locate  in  the  machine. 

Suitable  for  daily  use. 

5.  Fast  response: 

Response  time  for  producing  an  alarm  signal  should 
be  short  enough  to  prevent  the  following  teeth  being 
damaged . 

This  dissertation  deals  with  this  task  in  a 
deterministic  way  based  on  recognizing  and  analyzing  the 
tool  breakage  feature  in  a sensor  signal.  Then,  the  tool 
breakage  feature  is  extracted  by  a signal  processing 
algorithm  in  order  to  do  the  decision-making.  This  is 
different  from  the  stochastic  process  approach.  The 
contents  of  this  dissertation  are  summarized  below. 

In  Chapter  2,  the  basic  concepts  of  detection  of  tool 
breakage  throughout  the  work  are  discussed.  Considerable 
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research  has  been  undertaken  on  this  subject,  so  a 
literature  review  of  those  works  is  included. 

In  Chapter  3,  a mechanistic  model  of  the  milling 
process  is  described  and  computerized.  This  model  is  used 
to  simulate  the  dynamics  of  cutting  force  and  vibration  in 
milling  operations. 

In  Chapter  4,  a detection  algorithm  using  the  cutting 
force  signal  is  presented.  Experimental  work  to  verify 
this  algorithm  is  also  shown.  In  addition,  limitations  of 
cutting  force  signal  for  tool  breakage  detection  are 
discussed. 

In  Chapter  5,  tool  breakage  features  in  vibration 
signals  are  analyzed.  A detection  algorithm  with 
experimental  work  using  the  displacement  signal  is  shown. 

In  Chapter  6,  detection  of  tool  breakage  using  the 
displacement  signal  under  a variety  of  cutting  conditions 
is  shown. 

In  Chapter  7,  the  strategies  of  how  to  distinguish  the 
transient  response  signals  from  tool  breakage  signals  and 
how  to  set  an  alarm  threshold  are  explored  and  developed. 

Finally,  Chapter  8 concludes  this  work  and  presents 
some  ideas  about  future  research. 


CHAPTER  2 

BACKGROUND  AND  LITERATURE  REVIEW 


A system  of  detection  of  tool  failure  requires  very 
sensitive,  accurate,  and  reliable  techniques,  which  may  be 
classified  as  direct  and  indirect  methods. 

Direct  methods  generally  involve  measuring  the 
volumetric  change  in  the  tool  such  as  using  touch  trigger 
probes,  optical  sensors,  radio-active  sensors,  electrical 
resistance  sensors,  etc.  They  tend  to  be  off-line 
techniques  because  measurements  are  usually  taken  as  the 
tool  moves  out  of  the  cut.  Indirect  methods  utilize 
measurements  of  cutting  parameters  during  the  cutting  which 
may  be  correlated  with  the  tool  wear  or  breakage.  These 
cutting  parameters  include  temperature,  force,  torque, 
vibration,  sound,  acoustic  emission,  etc.  As  a result, 
indirect  methods  are  generally  considered  as  on-line 
techniques. 

As  mentioned  in  Chapter  1,  immediate  detection  of  tool 
failure  in  milling  should  be  implemented  in  real-time. 
Therefore,  direct  methods  may  not  be  appropriate  for  this 
application  because  they  tend  to  be  off-line  techniques. 
Optical  sensors  are  applicable  to  intermittent  cutting 
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operations  [5]  [6],  such  as  milling,  where  cutting  edges 
are  periodically  exposed  to  air.  Under  this  consideration, 
optical  sensors  may  be  utilized  as  an  on-line  technique  in 
milling.  However,  the  complexity  of  milling  operations  and 
workpiece  geometries  limits  its  application  as,  for 
example,  in  die  milling.  Further,  the  hostility  of  the 
cutting  environment,  such  as  dirt  and  lubrication  oil,  may 
lead  to  attenuation  and  scatter  of  light.  The  problem  of 
poor  reliability  and  repeatability  is  likely  to  be 
encountered  under  industrial  applications.  Therefore,  the 
detection  of  tool  breakage  in  milling  should  be  focused  on 
indirect  methods.  A general  basic  scheme  for  indirect 

methods  is  shown  in  Figure  2.1. 

Considerable  research  has  been  undertaken  using 
indirect  methods  [7].  In  recent  years,  two  approaches  have 
been  considered  to  be  worthy  of  further  investigation  and 
development.  One  approach  is  based  on  the  analysis  of  the 
acoustic  emission  signal  coinciding  with  tool  breakage; 
the  other  approach  is  based  on  the  analysis  of  changes  in 
the  pattern  of  the  cutting  force  signal. 


Figure  2 . 1 Tool  breakage  monitoring  scheme 
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2 . 1 Acoustic  Emission  Approach 

The  acoustic  emission  signal  is  obtained  from  a very 

high  frequency  piezoelectric  transducer.  Kannatey-Asibu  and 

Dornfeld  [8]  defined  the  acoustic  emission  as 

the  transient  elastic  energy  spontaneously 
released  in  materials  undergoing  deformation, 
fracture  or  both  [8,  pp.  248]. 

Therefore,  the  acoustic  emission  signal  is  sensitive  to 

micro-  and  macro-fractures  due  to  chip  formation,  chip 

breaking,  tool  wear,  and  tool  breakage  [9].  Typically,  the 

frequency  content  of  the  acoustic  emission  signal  is 

higher  than  100  KHz.  It  is  well  beyond  the  frequency  range 

of  machine  vibration  and  noise.  Thus,  the  relatively 

uncontaminated  acoustic  emission  signal  can  be  obtained  by 

passing  the  signal  through  a high-pass  filter.  However,  the 

choice  of  a suitable  bandwidth  of  the  acoustic  emission 

signal  varies  significantly  from  one  investigator  to  the 

next  [10]  [11]  [12]. 

It  has  been  shown  that  a sudden  acoustic  emission 
burst  is  generated  at  the  moment  of  tool  fracture  [13]. 
However,  the  acoustic  emission  signal  is  also  sensitive  to 
the  variation  of  cutting  parameters  [14].  Up  to  now,  no 
clear  relationship  has  been  established  between  the 
acoustic  emission  signal  and  the  cutting  parameters  [15]. 
Therefore,  it  is  not  easy  to  determine  the  proper  threshold 
value  to  distinguish  between  the  tool  breakage  and  the 
background  noise  due  to  chip  formation,  chip  breaking.  Two 
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recent  studies  have  been  proposed  to  solve  this  problem  in 
turning  where  single  point  machining  operations  are  used. 
One  is  a time  series  analysis  technique  using  an 
autoregressive  filter  to  model  the  acoustic  emission 
signal  [16].  Experiments  indicate  that  the  model 
parameters  could  ignore  the  change  of  cutting  parameters, 
but  show  a strong  sensitivity  to  the  progress  of  tool  wear. 
The  other  is  the  integration  of  cutting  force  and  acoustic 
emission  signals  via  neural  networks  [17].  The  cutting 
force  signal  was  measured  simultaneously  and  utilized  to 
solve  the  problem  of  threshold  setting  through  the  neural 
process  networks. 

However,  in  milling  operations,  a milling  cutter  with 
several  cutting  edges  can  lead  to  many  different  signal 
interactions.  The  ambiguous  identification  of  tool 
breakage  occurs  because  of  the  signal  interactions  in 
multi-tool  machining  [18].  In  addition,  entry  and  exit 
transients,  due  to  sudden  loading  and  unloading,  also 
produce  a large  acoustic  emission  burst  which  may  be 
higher  than  that  produced  by  tool  breakage,  depending  on 
the  extent  of  damage  on  the  tool  [19]  [20].  A special 
signal  processing  strategy  using  an  intelligent  insert  is 
being  developed  in  order  to  prevent  the  entry  and  exit 
spikes  from  interfering  with  tool  fracture  detection  [18]. 

To  summarize,  the  major  problems  of  the  acoustic 
emission  approach  appear  in  choosing  a suitable  bandwidth 
and  in  eliminating  the  noise  signals  due  to  chip 
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formation,  chip  breakage,  multi-tool  machining,  entry  and 
exit  transients,  the  variation  of  cutting  parameters,  etc., 
so  as  to  set  a proper  threshold  value.  In  addition,  the 
cost  of  the  required  instrumentation  is  also  expensive. 

2 . 2 Cutting  Force  Approach 

The  cutting  force  in  milling  offers  very  good  features 
for  sensing  the  milling  cutter  breakage.  There  are  various 
ways  of  measuring  the  cutting  force  in  milling,  either 
directly  from  a table-type  dynamometer,  or  a dynamometer 
built  into  the  spindle  bearings,  or  indirectly  evaluating 
the  cutting  force  from  the  spindle  motor  current,  torque, 
vibration,  or  sound.  Then,  the  cutting  force  signal  may  be 
processed  in  a number  of  ways  to  extract  features  depending 
on  the  desired  applications.  In  general,  the  processing 
algorithms  can  be  divided  into  either  the  time  series 
method  (black  box  type)  based  on  the  stochastic  process 
model  or  the  deterministic  method  based  on  the  analytical 
understanding  of  the  signal  changes. 

2.2.1  Time  Series  Method 

The  use  of  mathematical  models  to  explain  the 
behavior  of  physical  systems  is  the  foundation  of  science. 
If  the  mathematical  model  is  derived  by  the  physical  laws, 
the  model  is  called  a deterministic  model.  The  system 
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behavior  can  be  determined  or  calculated  accurately.  On  the 
other  hand,  if  the  mathematical  model  is  derived  by  the 
probabilistic  laws,  the  model  is  called  a stochastic  model. 
The  future  values  of  the  system  are  predicted  by  the 
stochastic  model  using  a sequence  of  time  series  data.  A 
general  process  can  be  described  by  the  combination  of  the 
two  models. 

In  a milling  process,  the  cutting  force,  f(t)  can  be 
represented  as  [20] 


Wg(B) 

0q(B) 

) = 

u(t)  + - 

a(t) 

(2-1) 

5r(B) 

♦ p(B) 

s = 

{Jo  “ W2^B  - cj2B^ 

**  • • • 

. - (JgB^ 

(2-2) 

r = 

1 - SiB  - 

• • • 

. - 

(2-3) 

q = 

1 - e^B  - 02B2 

• • • 

. - OgB^ 

(2-4) 

'P  = 

1 - <)>iB  - <^2®^ 

“ • • • 

. - ♦pBP 

(2-5) 

where 

B - backward  shift  operator 
u(t)  - cutting  process  parameters  which  are  assumed  to 
be  unknown 

a(t)  - gaussian  white  noise 

In  Equation  (2-1),  [(Js(B)/6i-(B)  ] u(t)  represents  the 
milling  dynamics  and  [eg(B)/«p(B) ] a(t)  represents  the 
process  noise.  In  order  to  make  the  detection  of  tool 
breakage  independent  of  cutting  process  parameters, 
[(js(B)/«Sr.(B)  ] u(t)  should  be  removed.  Thus  the  cutting 
force  versus  time  becomes  a stationary  stochastic  process; 
that  is. 
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0a(B) 

f(t)  = a(t)  (2-6) 

This  cutting  force  model  with  autoregressive  <^p(B)  and 
moving  average  9g(B)  components  is  called  an  autoregressive 
moving  average  model  (ARMA)  with  the  order  of  p and  q [21]. 
The  use  of  a real-time  recursive  ARMA(p,q)  model  requires 
nonlinear,  time  consuming  estimation  techniques.  However, 
an  autoregressive  model  (AR)  is  a linear  filter  and 
suitable  to  monitor  physical  systems  with  fast  dynamics. 
Therefore,  the  cutting  force  model  (AR)  with  the  order  p 
can  be  expressed  as 

f(t)  = a(t)  (2-7) 

<^p(B) 

= a(t)  (2-8) 

1 — <^2_B  — <^2B2  — • • • • — *^pBP 

and  then 

f(t)  = ♦if(t-l)  + «2f(t-2)  + + «^pf(t-p)  + a(t)  (2-9) 

The  value  of  f(t)  can  be  predicted  by  using  current 
estimated  values  of  ^2''  obtained  from  a 

recursive  parameter  adaptation  algorithm  and  past 
observation  f(t-l),  f(t-2),  that  is 

fi(t)  = «^i'f(t-l)  + <^2'f(t-2)  + ....  + <<>p'f(t-p)  (2-9) 

The  prediction  error,  a(t),  is  the  difference  between 

the  actual  f(t)  and  predicted  f ' (t) ; 

a(t)  = f(t)  - f(t)  (2-10) 

If  the  process  is  stationary,  the  prediction  error  is 
a white  noise  with  a constant  variance.  When  the  process 
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is  suddenly  deviated  by  a tool  breakage  event,  the  AR 
filter  is  disturbed.  A sudden  and  sharp  deviation  of 
prediction  error  is  observed.  This  is  how  the  time  series 
method  can  be  applied  to  tool  breakage  detection. 

In  [22],  a tool  breakage  detection  algorithm  has  been 
proposed  using  spindle  motor  current.  A 28th  order 
autoregressive  filter  was  applied  to  detect  tool  breakage. 
The  spindle  motor  current  is  known  to  be  proportional  to 
the  motor  torque  generated  in  the  cutting  process. 
Therefore,  the  changes  in  the  cutting  process,  for  instance 
tool  breakage,  can  be  monitored  by  using  the  spindle  motor 
current.  However,  experiments  demonstrated  that  the 
spindle  motor  current  is  proportional  to  the  cutting  torque 
only  up  to  about  20  Hz.  Thus,  this  algorithm  can  only  be 
applied  at  low  spindle  speeds. 

The  same  approach  was  also  implemented  using  feed 
force  in  end  milling  [23]  instead  of  spindle  motor  current. 
Although  the  cutter  had  only  two  flutes,  a 15th  order  of 
autoregressive  filter  was  still  required.  High  order  time 
series  filters  requiring  long  computation  time  windows 
might  not  be  adequate  for  real-time  applications.  Moreover, 
this  approach  still  could  not  distinguish  properly  between 
the  transient  cuttings  (i.e.  entry,  exit,  interrupted 
cuttings,  etc.)  and  tool  breakage  event. 
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2.2.2  Deterministic  Method 

It  seems  that  the  stochastic  process  approach  leads  to 
a very  high  order  time  series  modelling.  The  robustness  of 
detection  of  tool  breakage  is  still  not  satisfactory  with 
these  high  orders.  On  the  other  hand,  the  mechanics 
concerning  the  cutting  force  have  been  well  modelled  with 
regard  to  the  cutting  parameters  such  as  feedrate,  speed, 
depth  of  cut,  radial  immersion,  and  run-out.  Therefore,  the 
analytical  understanding  of  the  cutting  force  needs  to  be 
studied.  Once  the  basic  features  of  the  cutting  force  are 
recognized,  a more  robust  processing  algorithm  for 
detection  of  tool  breakage  is  possible. 

A scheme  for  monitoring  tool  and  cutting  conditions 
combining  five  features  extracted  from  the  cutting  force 
has  been  proposed  [24].  The  scheme  uses  these  five  basic 
features  and  is  able  to  separate  tool  breakage  from  other 
changes  in  the  cutting  process.  However,  the  scheme  works 
well  only  if  there  is  only  one  flute  actively  cutting.  With 
an  increasing  number  of  flutes  on  the  cutter,  it  becomes 
difficult  to  extract  these  five  features  from  the  cutting 
force. 

In  [25],  two  types  of  tool  breakage  detection  methods 
using  vibration  analysis  techniques  have  been  proposed.  One 
is  based  on  cepstrum  and  the  other  on  auto-correlation 
coefficients.  These  processing  techniques  require  a long 
finite-duration  time  window  with  a certain  time  delay.  For 
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the  transients  and  slotting,  a reliable  threshold  setting 
is  still  difficult. 

A tool  breakage  monitoring  method  using  the 
fluctuation  signal  of  the  spindle  rotational  speed  has  also 
been  proposed  [26].  The  fluctuation  signal  could  represent 
the  variation  of  the  cutting  torque  in  milling.  The 
detection  algorithm  is  based  on  a cutting  state  vector 
extracted  from  the  speed  fluctuation  pattern  during  one 
revolution.  Then,  the  cutting  state  vector  is  compared  with 
a reference  vector  to  recognize  tool  breakage.  The 
disadvantages  of  this  proposed  method  are  the  inherent  one 
revolution  delay  in  detecting  tool  breakage  and  the 
reference  vector  required. 

An  alternative  way,  using  the  average  force  per  tooth 
method,  has  been  developed  by  Altintas  et  al.  [27]  [28]. 
The  algorithm  is  based  on  synchronous  sampling  of  the 
cutting  force  with  the  spindle  speed  and  it  will  be  further 
developed  in  this  dissertation.  By  means  of  a simple 
pattern  analysis  [29],  tool  breakage  can  easily  be 
distinguished  from  transients  like  entry,  exit,  cornering, 
and  interrupted  cuttings. 

The  signal  used  as  the  basis  of  this  detection  method 
is  cutting  force  signal.  In  reality,  it  is  difficult  to 
accurately  measure  the  cutting  force  in  milling.  In  the 
experimental  work,  the  cutting  force  signal  is  measured 
using  a table-type  dynamometer,  but  it  distorts  the  cutting 
force  signal  at  high  frequencies  and  restricts  the  working 
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space  during  milling  operations.  Therefore,  compensation 
and  replacement  of  the  table-type  dynamometer  have  been 
presented  in  order  to  improve  these  problems  [30]. 
Alternative  signals,  like  an  accelerometer  signal  on  the 
headstock  and  a displacement  signal  of  the  spindle,  have 
been  considered  to  detect  tool  breakage  directly.  It  has 
been  shown  that  the  displacement  signal  correctly  indicates 
the  characteristic  changes  due  to  tool  breakage  [31].  Thus, 
a very  promising  way  for  the  detection  of  tool  breakage  in 
high  speed  milling  has  been  proposed  by  measuring  the 
displacement  of  the  spindle,  using  for  example  a 
capacitance  probe  [4]. 

In  addition,  a knowledge-based  approach  using  the 
average  displacement  signal  has  also  been  considered  [32]. 
It  is  expected  that  this  research,  adopting  the  knowledge- 
based  processing  with  AI  (artificial  Intelligence) 
technologies,  may  enhance  the  detection  system  which  has 
been  developed. 


CHAPTER  3 

MODELLING  OF  THE  MILLING  PROCESS 


3 . 1 Description  of  the  Milling  Process 

Milling  is  an  universal  machining  process  for  removing 
material  from  a workpiece  by  the  relative  motion  between 
the  workpiece  and  a rotating  cutter  having  multiple  cutting 
edges.  In  milling  operations,  the  surfaces  generating  by 
the  cutter  may  be  classified  into  two  general  categories, 
face  and  peripheral  millings. 

In  face  milling,  the  cutter  rotates  about  an  axis 
which  is  perpendicular  to  the  machined  surface  (see  Figure 

3.1.  a).  In  peripheral  milling,  the  cutter  rotates  about  an 
axis  which  is  parallel  to  the  machined  surface  (see  Figure 

3.1. b) . 

The  analysis  of  the  milling  process  has  been  well 
studied  [33-38].  Therefore,  a mechanistic  model  of  the 
milling  process  has  been  developed  [39-41]. 
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a) 


Figure  3.1  Types  of  milling  operations  [42,  pp.  148], 
a)  face  milling;  b)  peripheral  milling. 
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3 . 2 Milling  Force  Model 

In  a milling  operation,  the  tooth  path  resulting  from 
the  rotation  and  translation  of  the  cutter  is  a cycloid. 
However,  the  tooth  path  can  be  considered  as  circular  arcs 
with  good  approximation  under  the  assumption  that  the  feed 
per  tooth,  is  much  smaller  than  the  cutter  radius  (see 

Figure  3.2).  The  feed  per  tooth  is 

ft  = f / Nm  (3-1) 

where 

f - feedrate 
N - spindle  speed 

m - number  of  the  teeth  on  the  cutter 

In  Figure  3.2.b,  the  mean  chip  thickness,  h^,  at  the 
cutter  rotation  0 may  be  calculated  as 

hm  = ft  * sine  (3-2) 

3.2.1  Straight  Toothed  Cutter 

The  instantaneous  tangential  force  acting  on  the 
cutting  edge  is  the  product  of  the  specific  cutting  force 
and  the  area  of  the  chip  cut. 

fT  = ks  * 3a  * ^m 

= ks  * 3a  * ft  * sine  (3-3) 

where 

fij  - instantaneous  tangential  force  on  the 
cutting  edge 
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Figure  3.2  Geometry  of  milling  operation. 

a)  tooth  path;  b)  chip  thickness. 
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kg  - specific  cutting  force 
aa  - axial  depth  of  cut 

The  corresponding  instantaneous  radial  force  is  the 
instantaneous  tangential  force  multiplying  by  a constant, 

c . 

fj^  = c * f^  = c * kg  * a^  * ft  * sine  (3-4) 

To  decompose  the  cutting  forces  f-j  and  f^  into  the  X 
and  Y directions,  two  different  methods  of  generating 
surfaces  by  the  milling  operation  need  to  be  considered. 
They  may  be  called  up  milling  (see  Figure  3. 3. a)  and  down 
milling  (see  Figure  3.3.b). 

The  instantaneous  cutting  forces  in  the  X and  Y 
directions  can  be  represented  as 

fx  = fip  * COS0  + fj^  * sin0  (3-5) 

fY  = -fi  * sin0  + fR  * COS0  (3-6) 

In  multi-tooth  milling,  there  may  be  more  than  one 
tooth  active  in  the  cutting.  Therefore,  the  instantaneous 
cutting  forces  in  the  X and  Y directions  corresponding  to  a 
particular  angular  cutter  position  are  the  summation  of 
their  force  components  acting  on  each  tooth  engaged  with 
the  workpiece. 


m 

Fx  = S 

<S(i) 

* fx(i,0) 

(3-7) 

i=l 

m 

Fy  ~ 2 

5(i) 

* ^Y(i,0) 

(3-8) 

i=l 

and 
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Figure  3.3  Cutting  force  directions. 

a)  up  milling;  b)  down  milling. 


5(i)  = 1 
= 0 
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if  0s  < 0i  < 0Q 

otherwise 


where 

m - number  of  the  teeth  on  the  cutter 
0g  - starting  angle  of  cut 
0g  - exit  angle  of  cut 


3.2.2  Helical  Toothed  Cutter 


Figure  3. 4. a shows  a three  dimensional  helical 
toothed  cutter  and  the  incremental  cutting  force  along  the 
cutting  edges  of  the  cutter  with  the  axial  depth  of  cut, 
a^.  The  unrolled  circumference  of  the  cutter  is  shown  in 
Figure  3.4.b.  Every  helical  cutting  edge  becomes  a 
straight  line  under  the  helix  angle,  )3.  All  the  cutting 
edges  move  with  the  same  peripheral  cutting  speed.  The 
engagement  angle,  'If,  for  every  cutting  edge  is 

a^  * tan)S 

^ = (3-9) 

r 

Thus,  the  tooth,  T^,  moves  through  the  range  from  the 
position  0 = 0g  to  the  position  0 = 0^^  + v//.  Accordingly, 

the  teeth,  T2  and  T3 , move  by  'P  and  2'P  ahead,  respectively. 

Although  the  tool  may  have  the  several  teeth  cutting 

simultaneously,  the  actual  cutting  edges  exist  only  63  ®i 

< 0g . 

In  order  to  calculate  the  instantaneous  cutting 

force,  every  cutting  edge  is  considered  incrementally  by 


X 


Figure  3 


b) 


4 Geometry  of  a helical  toothed  cutter. 

a)  incremental  cutting  force; 

b)  unrolled  surface  of  cut. 
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, corresponding  to  the  increment  of  the  axial  depth  of 
cut,  Aa^.  Therefore,  the  element  of  tangential  force,  Af-p, 
is  the  product  of  the  specific  cutting  force  and  the 
corresponding  elemental  chip  cut. 

Afp,  = kg  * Aaa  * hj„ 


= kg  * Aa^  * ft  * sin0 
The  radial  component  of  force  is 

(3-10) 

Afj^  = c*Af.p  = c*kg*  Aa^  * f^  * sin0 

(3-11) 

Again,  these  elemental  forces  may  be  decomposed  into 
the  X and  Y directions.  The  elemental  forces  in  the  X and  Y 
directions  can  be  represented  as 


Afx  = Afip  * COS0  + Afj^  * sin0 

(3-12) 

Afy  = -Afpi  * sin0  + Afp  * COS0 

(3-13) 

To  obtain  the  instantaneous  force  acting 

on  the 

cutter  at  any  angular  position,  it  is  necessary  to 

sum  the 

elemental  forces  acting  along  the  engaged  ranges  of  all 
the  cutting  edges. 


m n 

Fx  = 2 S 6(i,j)  * Afx(i,j) 

i=l  j=l 

(3-14) 

m n 

Fy  = Z Z (S(i,j)  * AfY(i^j) 

i=l  j=l 

and 

5(i,j)  =1  if  0g  < 6(i,j)  < 9e 
= 0 otherwise 

(3-15) 

where 

n - total  incremental  elements  in  the  direction  of 


the  axial  depth  of  cut 
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m - number  of  the  teeth  on  the  cutter 

In  the  foregoing  discussion,  it  has  been  assumed  that 
the  cutting  force  is  obtained  from  a cutter  designed  with 
equally  spaced  cutting  edges,  with  no  consideration  of  the 
regeneration  of  waviness  and  the  run-out  on  the  cutter.  In 
reality,  both  of  them  not  only  exist  but  also  vary  the 
cutting  force.  Therefore,  they  should  be  taken  into  account 
in  the  cutting  force  model. 

3 . 3 Regeneration  of  Waviness 

Because  of  the  inherent  insufficient  stiffness  of  the 
spindle  system,  there  is  a cutter  vibration  at  the  time 
that  the  workpiece  surface  is  produced.  As  a result,  a 
waviness  is  generated  on  the  workpiece  surface  [38]  (see 
Figure  3. 5. a).  A variable  component  of  chip  thickness,  h^, 
is  obtained  due  to  the  regeneration  of  the  waviness.  This 
is  because  each  tooth  removes  the  chip  from  the  surface 
produced  by  the  preceding  teeth.  If  any  vibration  exists 
during  the  cut,  a variable  thickness  is  generated. 

hv  = z - Zb  (3-16) 

where 

z - current  vibration  of  the  cutter  normal  to  the  cut 
Zb  - undulation  of  the  surface  produced  by  the 
preceding  teeth  at  current  position  of  cut 
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a) 


Figure  3.5  Regeneration  of  the  waviness  in  milling. 

a)  waviness  of  workpiece  surface; 

b)  simulation  model. 
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A milling  cutter  model  may  be  established  by  using 
modal  analysis  on  the  transfer  function  of  the  spindle 
system  (see  Figure  3.5.b).  The  vibration  of  the  cutting 
edge  may  be  expressed  as 

z = X * sin0  + y * cos0  (3-17) 

where 

X - displacement  of  the  cutter  in  the  X direction 
y - displacement  of  the  cutter  in  the  Y direction 
Once  the  instantaneous  cutting  force  has  been 
calculated,  the  vibrations  of  x^  and  y^  in  each  mode  can 
also  be  solved  by  Euler's  equations. 


“ ^Xi  * - Kxi  * Xi 

Mxi 


(3-18) 


Yj  = 


Fy  - Cyj  * yj  - Kyj  * yj 

Myi 


(3-19) 


and 


Xi(t+At)  = Xi(t) 

+ Xi(t)  * At 

(3-20) 

yj (t+At)  = yj (t) 

+ yj (t)  * At 

(3-21) 

Xi(t+At)  = Xi(t) 

+ Xi(t+At)  * At 

(3-22) 

yj (t+At)  = yj (t) 

• + Yj (t+At)  * At 

(3-23) 

where 

x^,  yj  - acceleration  of  the  X or  Y direction  in  i 
or  j mode 

Xi,  yj  - velocity  of  the  X or  Y direction  in  i or 


j mode 
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Xi,  Yj  - displacement  of  the  X or  Y direction  in  i 
or  j mode 

Mxi,  Myj  - mass  coefficients 
Cxi/  ^Yj  “ damping  coefficients 
Kxi,  Kyj  - spring  coefficients 
At  - integration  time  interval 

Finally,  the  vibrations  of  cutter  in  the  X and  Y 
directions  can  be  obtained  by  the  superposition  of  the 

vibration  in  each  mode. 

nx 

X = S Xi  (3-24) 

i=l 

ny 

y = S Yj  (3-25) 

j=l 

where 

nx  - number  of  modes  in  the  X direction 

ny  - number  of  modes  in  the  Y direction 

Substituting  Equations  (3-24)  and  (3-25)  into  (3-17) , 
the  actual  chip  thickness,  h,  can  be  expressed  as  the 

summation  of  the  mean  and  variable  components  of  chip 
thickness.  That  is 

h = hjn  + hv  = ft  * + z - Zb  (3-26) 

This  is  only  true  when  (ft  * sin0  + z)  > Zb*  If  the 
vibration  between  the  cutter  and  workpiece  becomes  so  big 
that  the  cutting  edge  may  jump  out  and  loose  contact  with 
the  workpiece  surface,  i.e.  (f^  * sin0  + z)  < z-^,  the 

actual  chip  thickness,  h,  is  equal  to  zero. 
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Now,  Equations 

(3-3) 

and  (3-10)  may 

be  modified  to 

incorporate  the  effect  of 

the  regeneration 

of  waviness  and 

become 

fT  = ks  * 

aa  * 

(ft  * sin0  + z - 

Zb) 

if 

ft  * 

sin0  + z > zjj 

= 0 

if 

ft  * 

sin0  + z < zj3 

(3-27) 

and 

fT  = ks  * 

Aaa 

* (ft  * sin0  + z 

- Zb) 

if 

ft  * 

sin0  + z > zt 

= 0 

if 

ft  * 

sin0  + z < Zj3 

(3-28) 

Finally,  it  has 

been 

shown  that  the 

influence  of  the 

regeneration  of  waviness  on  the  cutting  force  may  be 
neglected  as  long  as  milling  is  in  the  stable  cut  without 
chattering  [43]. 


3.4  Run-out  on  the  Cutter 

Cutter  run-out  is  a common  phenomenon  affecting  the 
performance  of  milling  operations.  In  general,  run-out  may 
be  referred  to  as  radial  and  axial  throws  which  may  arise 
from  eccentricities  when  the  cutter  is  held  in  the  spindle. 
When  indexable  insert  type  cutters  are  used,  run-out  is 
more  likely  due  to  insert  size  variation,  spindle  out-of- 
roundness,  and  setting  of  inserts  [44].  For  the  non- 
indexable,  regrindable  type  cutters,  run-out  commonly 
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results  from  the  cutter  offset  with  respect  to  the  spindle 


rotation  [41]. 

Usually,  the  axial  throw  is  small  relative  to  the 
axial  depth  of  cut,  so  its  effect  on  the  cutting  force  may 
be  neglected.  However,  the  magnitude  of  radial  throw  is 
the  same  order  as  that  of  feed  per  tooth.  It  may  cause 
each  tooth  to  experience  a different  chip  thickness. 
Therefore,  the  cutting  force  analysis  becomes  much  more 
complicated  due  to  the  presence  of  the  radial  throw. 

The  chip  thickness  with  radial  throw  can  be  derived  by 
using  the  circular  tooth  path  approximation  [41],  see 
Figure  3.6.  The  actual  mean  chip  thickness,  hj^^f  can  be 
expressed  as 


^mc  “ ^i  ” ^b 


(3-29) 


where 


rj^  - radius  of  the  ith  tooth  in  the  cut 
rjj  - distance  from  the  ith  tooth  center  to  surface 
left  by  the  preceding  cut 
Using  the  law  of  sines,  r^,  can  be  calculated  as 

^i-i 


sin(7T/2  + 0) 


sin(7T/2  - 0) 


(3-30) 


then 


COS0 

= ^i-1  

COS0 


(3-31) 


and 


0 = 0 + e 


(3-32) 
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Figure  3.6  Geometry  of  chip  thickness  with  radial  throw. 
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Substituting  Equation  (3-32)  into  (3-31) , 

sin0  * sine 

rjj  = * (cose ) (3-33) 

COS0 


Again,  using  the  law  of  sines,  cose  and  sine  can  be 
calculated  as 


sin(7r/2  + 0)  sine 

then 


and 


sine 


* COS0 


^i-1 


(3-35) 


(ri_l  - * cos20)i/2 

cose  = (1  - sin^e)^/2  = (3-36) 

^i-1 

Substituting  Equations  (3-35)  and  (3-36)  into  Equation 
(3-33) , can  be  expressed  as 

rj3  = * cos^0)V2  - * sin0  (3-37) 

Finally,  substituting  Equation  (3-37)  into  (3-29),  hmc 
can  be  expressed  as 

hjnc  = ~ ^t^  * cos20)V2  + * sin0  (3-38) 

The  term,  (ft^  * cos20),  may  be  neglected  in  Equation  (3- 
38) , since  the  feed  per  tooth,  f^,  is  much  smaller  than  the 
cutter  radius,  i.e.  >>  ft^  * cos20)  . Therefore,  the 

actual  mean  chip  thickness  is 

hmc  = (^i  " ^i-l)  + ft  * si>30  (3-39) 

This  result  shows  that  radial  throw  introduces 
another  term  to  the  chip  thickness.  It  also  needs  to  be 
considered  that  the  radial  throw  may  be  severe  in 
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comparison  to  the  term,  ( * sin0),  especially  at  entry 
and  exit  positions  where  0 is  small.  As  a result,  r^.j^  may 
be  greater  than  (r^  + * sin0)  which  results  in  negative 

chip  thickness  value  in  Equation  (3-39) . What  is  really 
happening  is  that  some  teeth  experience  no  chip  load  or 
partial  chip  load  and  a particular  tooth  removes  the 
workpiece  surface  left  two  or  more  teeth  ago. 

Therefore,  Equation  (3-39)  needs  a further 
modification  to  include  these  conditions. 

^mc  = (^i  " ^i-k)  + k * ft  * sin0  (3-40) 

where  tooth  i removes  the  surface  left  by  tooth  i-k. 

However,  the  determination  of  the  coefficient,  k,  is 
not  straightforward.  Hence,  a special  calculation  algorithm 
has  been  developed  [39].  It  checks  the  previous  chip 
thickness  consecutively  until  the  constant,  k,  is  found. 
That  is 

k 

hi  = S ft  * sin0  + ri_j  - 

j=l 

if  ft  * sin0  + ri_]^  - < 0 (3-41) 

where  h^  is  a non-positive  constant,  h^  may  also  be  called 
the  chip  thickness  of  over  cut  due  to  the  radial  throw. 

Hence,  the  modified  actual  mean  chip  thickness  may  be 
expressed  as 

hmc  = ft  * sin0  + (ri  - ri_i)  + hi 
if  ft  * sin0  + ri  > ri_i 
= 0 


if  ft  * sin0  + ri  < ri_i 


(3-42) 
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Again,  Equations  (3-3)  and  (3-10)  may  be  modified  to 
incorporate  the  effect  of  radial  throw  and  become 
fT  = kg  * aa  * (ft  * sine  + ri  - ri_i  + hi) 
if  ft  * sine  + ri  > ri_i 

f iji  = 0 

if  ft  * sine  + ri  < ri_i  (3-43) 

and 

ftp  = kg  * Aa^  * (ft  * sine  + ri  - ri_p  + hi) 
if  ft  * sine  + ri  > ri_p 

f iji  = 0 

if  ft  * sine  + ri  > ri_p  (3-44) 

Finally,  Figure  3.7  shows  a flow  chart  of  modelling 
the  dynamics  of  cutting  force  in  milling.  A mechanistic 
force  model  with  straight  and  helical  toothed  cutters 
including  the  effects  of  regeneration  of  waviness  and 
radial  throw  has  been  computerized,  see  Appendix  [A]  . In 
the  following  chapters,  this  model  will  be  utilized  as  a 
useful  tool  for  the  simulation  and  analysis  of  the  cutting 
force  signal  in  milling  operations. 


Figure  3.7  Flow  chart  for  simulation  of  cutting  force 


CHAPTER  4 

DETECTION  OF  TOOL  BREAKAGE  USING  THE  CUTTING  P'ORCE  SIGNAL 


4 . 1 Simulation  of  the  Cutting  Force  Signal 
With  and  Without  Tool  Breakage 

4.1.1  Straight  Toothed  Cutter 

The  cutting  force  signal  in  face  milling  is  periodic 
with  tooth  frequency,  ff  As  one  of  the  teeth  is  chipped 
or  broken,  the  tooth  following  the  chipped  or  broken  one 
will  have  a bigger  chip  load  in  order  to  remove  the 
material  left  by  the  chipped  or  broken  one.  Therefore,  the 
characteristics  of  the  cutting  force  signal  are  changed 
because  of  the  existence  of  tool  breakage. 

Figure  4.1  shows  a face  milling  cutter  (diameter  = 
101.6  mm)  without  run-out  milling  over  the  most  common 
geometry  of  workpiece.  The  cutting  parameters  are  shown  in 
Table  4.1.  Based  on  the  transfer  functions  measured  [45], 
the  model  of  the  spindle  system  at  the  end  of  tool  is  shown 
in  Table  4.2.  Figure  4 . 2 shows  the  simulation  of  the 
cutting  force  signal  in  the  entry,  steady,  and  exit  states 
with  and  without  tool  breakage.  The  simulation  results  show 
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Figure  4.1  Cut  geometry  in  the  entry,  steady, 
and  exit  states. 
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Undamaged  cutter 

2.0 
1.0 
0.0 


Damaged  Cutter 


Figure  4.2  Simulation  of  the  cutting  force  in  the  entry, 
steady,  and  exit  states  in  up  milling. 
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the  effect  of  regeneration  of  waviness  on  the  cutting  force 
is  slight  and  can  be  neglected  because  the  cut  is  stable. 

By  experiments,  it  was  seen  that  the  broken  tooth 
breaking  off  the  whole  cutting  area  was  common.  Therefore, 
it  was  assumed  that  the  cutting  area  of  a broken  tooth  was 
totally  broken  in  the  simulation. 

Table  4.1  Cutting  parameters  in  the  simulation  of 
cutting  force  in  face  milling 


spindle 

depth  of 

feed  per 

radial 

workpiece 

speed 

cut 

tooth 

immersion 

material 

rpm 

mm 

mm 

600 

1.27 

0.254 

1/2 

cast  iron 

Table  4.2  Simulation  parameters  for  the  spindle  system 


direction 

mass 

damping 

natural 

(kg) 

ratio 

frequency 

X 

18.6 

0.04 

468  Hz 

Y 

19.2 

0.04 

515  Hz 

In  the  entry  transient,  the  cutter  gradually  enters 
the  workpiece  and  the  width  of  the  cutting  force  signal 
increases.  As  to  the  steady  state,  the  cutting  force  signal 
is  periodic  per  tooth  period  because  of  the  invariant 
cutting  geometry.  In  the  exit  transient,  both  the  magnitude 
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and  the  width  of  the  cutting  force  signal  decrease. 
However,  for  the  damaged  cutter  with  a broken  tooth,  the 
cutting  force  signal  which  decreases,  then  increases,  forms 
a special  tool  breakage  feature  shown  in  all  of  the  states. 
It  is  also  true  for  the  case  of  down  milling,  see  Figure 
4.3. 

Analogous  results  were  also  obtained  from  the  other 
immersions.  Figure  4.4  shows  the  cutting  force  signal  in 
the  steady  state  with  1/8,  1/4,  1/2,  and  1/1  immersions  of 
cut  (where  1/1  immersion  represents  slotting) . In  addition. 
Figure  4.5  shows  the  cutting  force  signal  in  the  steady 
state  using  the  cutter  with  4,  6,  8,  and  10  teeth  in  the 
1/2  immersion.  The  clear  tool  breakage  feature  is  available 
in  all  of  the  cases. 

Run-out  is  an  unavoidable  phenomenon  in  milling.  In 
experiments,  the  best  radial  throw  on  the  cutting  edges  was 
found  to  be  within  0.025  mm  (0.001”).  An  example  of  radial 
throw  on  each  tooth  is  listed  in  Table  4.3. 

As  mentioned  in  Chapter  3 , the  existence  of  run-out 
produces  the  variation  of  chip  load  for  each  tooth. 
Therefore,  the  variable  chip  load  affects  the 
characteristics  of  the  cutting  force  signal.  A simulation 
of  the  cutting  force  signal  in  entry,  steady,  and  exit 
states  with  these  radial  throws  is  shown  in  Figure  4.6.  The 
cutting  parameters  are  the  same  as  those  in  Figure  4.2.  It 
can  be  seen  that  the  periodicity  of  the  cutting  force 
signal  per  tooth  is  disturbed  by  the  radial  throws. 
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Undamaged  Cutter 


Damaged  Cutter 


Figure  4.3  Simulation  of  the  cutting  force  in  the  entry, 
steady,  and  exit  states  in  down  milling. 
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Undamaged  Cutter 
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2.0 


1.0 
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Figure  4.4  Simulation  of  the  cutting  force  in  the  1/8, 
1/4,  1/2,  and  1/1  immersions  of  cut. 
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Undamaged  Cutter  Damaged  cutter 


Figure  4.5  Simulation  of  the  cutting  force  with  4,  6,  8, 
and  10  teeth  on  the  cutters. 
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Undamaged  cutter 
* 1o3(n)  a) 


Damaged  Cutter 


* lO^(N)  d) 


Figure  4.6  Simulation  of  the  cutting  force  with  radial 
throw  in  the  entry,  steady,  and  exit  states. 
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However,  much  bigger  variation  due  to  tool  breakage  still 
remains. 


Table  4.3  Radial  throw  in  the  face  milling  cutter 


tooth 

number 

radial  throw 
(mm) 

tooth 

number 

radial  throw 
(mm) 

1 

-0.0203 

5 

-0.0252 

2 

-0.0042 

6 

-0.0029 

3 

-0.0160 

7 

0.0 

4 

-0.0080 

8 

-0.0069 

4.1.2  Helical  Toothed  Cutter 


In  peripheral  milling,  one  of  the  most  common  types  is 
end  milling.  Figure  4.7  shows  the  simulation  of  the  cutting 
force  signal  without  run-out  in  the  1/2  immersion  up 
milling  by  using  a 4-flute  end  mill  cutter  (diameter  = 
25.4  mm)  with  a helix  angle  of  30°.  The  cutting  parameters 
are  listed  in  Table  4.4. 

Even  though  the  cutting  edges  are  helical,  the 

I 

variation  of  the  cutting  force  signal  is  still  periodic 
with  the  tooth  frequency.  The  uniform  periodicity  is 
disturbed  when  one  of  the  cutting  edges  is  broken.  As  a 
result,  much  bigger  variation  of  the  cutting  force  signal. 


46 


Und2Uttaged  Cutter 
* 1o3(N) 


b) 


Damaged  Cutter 

* 1o3(N)  d) 

4.O1 , 


f) 


Figure  4.7  Simulation  of  the  cutting  force  in  the  entry, 
steady,  and  exit  states  using  a helical 
toothed  cutter. 
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which  decreases,  then  increases,  is  generated  because  of 
the  broken  edge. 


Table  4.4  Cutting  parameters  in  the  simulation  of  the 
cutting  force  in  end  milling 


spindle 

depth  of 

feed  per 

radial 

workpiece 

speed 

cut 

tooth 

immersion 

material 

rpm 

mm 

mm 

600 

2.54 

0.254 

1/2 

cast  iron 

This  is  also  true  for  the  different  immersions  of 
cut.  Figure  4.8  shows  the  cutting  force  signal  in  the  1/8, 
1/4,  1/2,  and  1/1  immersions  of  cut.  Figure  4.9  shows  the 
cutting  force  signal  with  2,  3,  and  4 flutes  of  the  end 
mill  cutter  in  the  steady  state.  It  is  shown  that  the 
distinction  between  the  undamaged  and  damaged  cutters  is 
stronger  when  a less  number  of  flutes  on  the  cutter  is 
used. 

In  the  foregoing  discussion,  it  is  clearly  seen  that 
the  tool  breakage  feature  in  the  cutting  force  signal  is 
available  in  all  kinds  of  milling  operations.  Therefore, 
processing  the  cutting  force  signal  is  recognized  as  one  of 
the  most  feasible  way  to  detect  tool  breakage  in  milling. 
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Undamaged  cutter 

4.0 

2.0 
0.0 


* 1o3(N) 
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Damaged  Cutter 


Figure  4.8  Simulation  of  the  cutting  force  in  the  1/8, 
1/4,  1/2,  and  1/1  immersions  of  cut  using  a 
helical  toothed  cutter. 
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Undamaged  Cutter 


Damaged  Cutter 


(SEC) 


Figure  4.9  Simulation  of  the  cutting  force  using 
helical  toothed  cutters  with  2,  3, 
and  4 flutes . 
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4 . 2 Detection  Algorithm  for  Tool  Breakage 

In  Section  4.1,  it  has  been  shown  that  the 
characteristics  of  the  cutting  force  signal  are  obviously 
changed  due  to  tool  breakage.  Hence,  the  cutting  force 
signal  is  available  for  recognizing  tool  breakage.  However, 
it  is  difficult  to  detect  tool  breakage  using  the 
instantaneous  cutting  force  signal  directly.  This  is 
because  cut  geometry,  run-out,  and  the  difficulties  of 
force  measurement  (Section  4.5)  change  and  distort  the 
cutting  force  signal.  Therefore,  it  is  necessary  to  develop 
a detection  algorithm,  and  extract  the  tool  breakage 
feature  from  the  instantaneous  cutting  force  signal. 

In  order  to  understand  the  tool  breakage  feature  in 
the  cutting  force  signal,  the  spectrum  of  the  cutting  force 
signal  is  utilized.  Figure  4.10  shows  the  cutting  force 
spectra  corresponding  to  the  signal  in  Figure  4.2.  The  DC 
component  of  the  spectrum  is  omitted  because  the  tool 
breakage  feature  belongs  to  the  AC  components  in  the 
cutting  force  spectrum.  For  the  undamaged  cutter,  the 
prominent  components  are  not  only  the  component  of  tooth 
frequency,  but  also  its  harmonics.  Both  are  needed  to 
represent  the  periodic  change  of  the  cutting  force  signal 
exactly.  For  the  damaged  cutter  with  a broken  tooth,  the 
component  of  once  per  revolution,  fj-,  and  its  harmonics 
show  up.  In  comparing  the  spectra  between  the  undamaged  and 
damaged  cutters,  it  is  clearly  shown  that  the  frequency 
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Figure  4.10  Spectra  of  the  simulation  of  the  cutting 

force  in  the  entry,  steady,  and  exit  states. 
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components  of  tool  breakage  are  and  its  harmonics, 

depending  on  the  number  of  teeth  on  the  cutter.  This 
frequency  area  may  be  called  a "tool  breakage  zone" . 

Now,  the  task  is  how  to  process  the  cutting  force 
signal  and  make  the  components  inside  the  tool  breakage 
zone  stand  out,  or  suppress  the  components  outside  this 
zone.  Several  methods  have  been  under  investigation,  such 
as  the  time  series  method,  the  power  method  [46],  and  the 
difference  per  tooth  average  method.  So  far,  it  seems  that 
the  difference  per  tooth  average  method  is  the  best  one. 
The  main  advantage  of  this  method  is  that  it  is  simple, 
quick,  and  reliable.  The  difference  per  tooth  average 
method  is  illustrated  in  the  following  text: 

The  average  force  per  tooth  is  calculated  as 

n 

S Fx(i) 
i=l 

FxavgC:)  - 


yavg 


(j) 


n 


S Fy(i) 
i=l 


n 


(4.2) 


where 

Fjjavg  “ average  force  per  tooth  in  the  X direction 
Fyavg  “ average  force  per  tooth  in  the  Y direction 
j - index  for  time  increment  per  tooth  period 
n - number  of  samples  per  tooth  period 
In  reality,  machining  operations  may  involve  contours 
and  corners.  Hence,  both  the  X and  Y directions  of  the 
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cutting  force  should  be  taken  into  account  and  calculated 
as  the  resultant  average  force  per  tooth,  F^vg- 

^avg(^)  “ (^xavg^(j)  "*■  ^yavg^  ( ^ (4-3) 

Finally,  the  difference  of  the  resultant  average  force 
per  tooth,  AF^^g  is  obtained  by 

^^avg(^)  ~ ^avg(j)  ~ ^avg(j“^)  (4.4) 

Applying  the  difference  per  tooth  average  method. 
Figure  4.11  shows  the  simulation  results  for  the  cutting 
force  signal  in  the  steady  state  as  shown  in  Figure  4.2. 
The  first  row  displays  the  cutting  force  signal.  For 

calculation  of  the  average  force,  nine  samples  were  taken 
in  each  tooth  period.  Averaging  the  cutting  force  signal 
per  tooth,  see  Equations  4.1  and  4.2,  filters  out  the  tooth 
frequency  and  its  harmonic  components.  Then,  the  resultant 
average  cutting  force  signal  per  tooth,  see  Equation  4.3, 
is  displayed  in  the  second  row.  It  is  shown  that  the 

periodic  variation  per  tooth  is  removed.  However,  the  DC 
component  of  the  cutting  force  and  the  variation  due  to 
the  broken  tooth  still  remain.  The  variation  due  to  the 
broken  tooth  stands  out,  once  the  average  of  the  cutting 
force  per  tooth  is  differenced,  see  Equation  4.4.  The 

difference  of  the  average  cutting  force  signal  is 
displayed  in  the  third  row.  It  is  shown  that  there  is  an 
extremely  clear  distinction  of  the  difference  of  the 

average  cutting  force  between  the  undamaged  and  damaged 
cutters.  Therefore,  this  processing  method  is  a very 
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Figure  4.11  Simulation  of  the  instantaneous  force, 
average  force,  and  difference  of  the 
average  force  in  the  steady  state. 
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efficient  way  of  extracting  the  tool  breakage  feature  from 
the  cutting  force  signal. 

This  is  true  not  only  in  the  steady  state,  but  also  in 
the  transient  states.  Figure  4.12  and  Figure  4.13  show  the 
simulation  results  in  the  entry  and  exit  transients.  A 
staircase  of  the  average  force  signal  is  obtained  due  to 
the  finite  sampling.  However,  the  staircase  will  disappear 
when  a higher  sampling  rate  is  used  [27]  [29]. 

Run-out  in  the  form  of  radial  throw  also  produces 
another  kind  of  once  per  revolution  signal,  see  Figure 
4.14.  Figure  4.14.b  and  4.14.f  show  the  spectra  of  the 
cutting  force  data  with  radial  throw  in  the  steady  state 
portion  of  cut  indicated  in  Figure  4.6.  Unfortunately,  the 
frequency  components  of  the  radial  throw  are  also  located 
inside  the  tool  breakage  zone.  Therefore,  it  affects  the 
difference  of  the  average  force.  Figure  4.14.C  shows  the 
average  force  signal  with  the  radial  throw.  Even  for  the 
undamaged  cutter,  a certain  level  of  the  difference  of  the 
average  cutting  force  signal  exists  (see  Figure  4.14.d). 
However,  a clear  distinction  of  the  difference  of  the 
average  signals  is  still  available  between  the  undamaged 
and  damaged  cutters. 

To  summarize,  processing  the  cutting  force  signal 
with  the  difference  per  tooth  average  method  through 
simulation  is  recognized  as  a very  feasible  way  to  detect 
tool  breakage  in  milling.  Thus,  an  experimental 
verification  of  this  method  is  necessary. 
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Figure  4.12  Simulation  of  the  instantaneous  force, 
average  force,  and  difference  of  the 
average  force  in  the  entry  state. 
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Figure  4.13  Simulation  of  the  instantaneous  force, 
average  force,  and  difference  of  the 
average  force  in  the  exit  state. 
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Figure  4.14  Simulation  of  the  instantaneous  force  and  its 
spectrum,  average  force,  and  difference  of 
the  average  force  in  the  steady  state  with 
radial  throw. 
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4 . 3 Experimental  Set-up 

A Sundstrand  Omnimil  20  machining  center  (White 
Sundstrand  Machine  Tool  Co.)  shown  in  Figure  4.15,  with  a 
maximum  rated  spindle  speed  of  5500  rpm  and  a rated  horse 
power  of  25  at  1200  rpm  was  used  in  the  experiments.  It  is 
a horizontal  spindle  machine  with  three  linear  axes, a 360 
discrete  degree  rotary  index  table,  and  a 30  position 
automatic  tool  exchanger  [46]. 

All  of  the  face  milling  tests  were  done  by  using  a 4- 
inch  diameter  double-negative  ZN3M  Carbaloy  milling  cutter 
(GE  Co.).  The  catalog  number  of  the  inserts  is  SNG634  or 
SNG632  (silicon  nitride) . The  maximum  run-out  among  these 
inserts  on  the  cutter  was  0.0254  mm  (0.001").  The  workpiece 
material  was  an  ingot  of  GM-241  cast  iron  in  the  form  of  a 
flat  12"  long,  6”  wide  and  3.5”  thick.  The  workpiece  was 
bolted  to  an  angle  plate  so  that  a large  surface  was 
available  for  milling. 

The  cutting  force  signal  was  obtained  by  using  a 
table  type  of  dynamometer  under  the  workpiece  and  the  angle 
plate,  see  Figure  4.16.  The  dynamometer  was  constructed 
with  three  directional  Kistler  load  cells  (Model  9067)  . The 
output  signals  of  the  load  cells  travel  through  wires  , to 
Kistler  charge  amplifiers  (Type  5007) . The  charge 
amplifiers  convert  the  signals  yielded  by  the  load  cells 
into  the  proportional  cutting  force  signals  in  three 
perpendicular  directions.  The  transfer  function  of  the 
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Figure  4.15  Series  20  Omnimil  Machining  Center. 
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Figure  4.16  Assembly  of  workpiece-angle  plate-dynamometer 
and  transfer  function  of  the  dynamometer, 
a)  Y direction;  b)  Z direction. 
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dynamometer  in  the  Y and  Z directions  is  shown  in  Figure 
4. 16. a and  4.16.b.  This  transfer  function  indicates  the 
accuracy  with  which  the  dynamometer  signal  represents  the 
actual  cutting  force  signal  (see  Section  4.5).  The  cutting 
force  signals  were  recorded  on  an  IBM  PC-AT  microcomputer 
using  a data  acquisition  board  (Data  Translation  DT2818). 
The  data  acquisition  system  provides  four  A/D  converters 
with  a maximum  simultaneous  sampling  rate  of  27.5  KHz  and 
two  D/A  converters  with  a maximum  sampling  rate  of  33.0 
KHz.  The  voltage  resolution  of  this  board  is  +/-  10  VDC 
with  12  bits. 

In  order  to  sample  the  cutting  force  signal 
synchronized  with  the  spindle  rotation,  external  sampling 
and  triggering  were  used.  A magnetic  pick-up  sensed  the 
distance  between  itself  and  a steel  gear  with  72  teeth 
mounted  on  the  spindle.  When  the  spindle  rotated,  the 
distance  was  converted  into  an  analog  signal.  This  analog 
signal  was  transmitted  through  a comparator,  a Schmitt 
trigger,  and  a counter  from  which  the  external  sampling  and 
triggering  signals  were  obtained.  Figure  4. 17. a shows  the 
electrical  circuit  block  diagram  for  generating  the 
external  sampling  and  triggering  signals.  Figure  4.17.b 
shows  the  interrupt  circuit  for  feed  control  of  the 
machining  center  when  tool  breakage  is  detected.  Finally 
the  schematic  diagram  for  the  experimental  work  is  shown  in 
Figure  4.18. 
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gure  4.17  Electrical  circuit  block  diagram. 

a)  external  sampling  and  triggering; 

b)  interrupted  circuit  for  feed  control. 
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Figure  4 . 18  Process  diagram  for  the  experimental  work 
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4 . 4 Experimental  Verification 

Using  the  experimental  set-up  described  in  Section 
4.3,  a series  of  cutting  tests  was  carried  out.  They 
correspond  to  the  simulation  results  shown  in  Sections  4.1 
and  4.2. 

Figure  4.19  shows  a 1/2  immersion  up  milling  case.  The 
cutting  parameters  are  the  same  as  listed  in  Table  4.1.  A 
insert  with  a broken  cutting  edge  was  used  to  obtain  the 
tool  breakage  signal  in  the  whole  cycle  of  cutting  tests. 
The  first  row  displays  the  instantaneous  cutting  force 
signal  in  the  whole  cycle.  The  average  force  signal  is 
displayed  in  the  second  row  and  the  difference  of  the 
average  force  signal,  in  the  third  row.  A clear 
distinction  is  available  between  the  undamaged  and  damaged 
cutters.  In  the  transient  states,  the  staircase  shape  of 
the  simulation  of  the  average  force  disappears.  It  is 
smoothed  not  only  by  the  existence  of  run-out  but  also 
further  by  the  characteristics  of  the  transfer  function  of 
the  dynamometer.  Figure  4.20  is  an  expanded  view  of  Figure 
4.19  during  the  steady  state  portion  of  the  cut.  The  first 
row  displays  the  instantaneous  cutting  force  signal.  The 
second  row  displays  the  corresponding  spectra  of  the 
cutting  force  signal.  The  third  row  displays  the  average 
force  signal;  the  fourth  row,  the  difference  of  the  average 
force  signal.  The  feature  distinguishing  tool  breakage  is 
clearly  observed  in  the  time  domain  and  also  in  frequency 
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Figure  4.19  Instantaneous  force,  average  force,  and 
difference  of  the  average  force  in  the 
complete  cutting  cycle. 
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Figure  4.20  Instantaneous  force  and  its  spectrum,  average 
force,  difference  of  the  average  force  in  the 
steady  state. 
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domain  of  the  tool  breakage  zone.  The  average  force  and  the 
difference  of  the  average  force  signal  for  the  1/8,  1/4, 
3/4  and  1/1  are  also  shown  in  Figure  4.21  and  Figure  4.22. 

By  experiments,  it  is  shown  that  applying  the 
difference  per  tooth  average  method  to  process  the  milling 
force  signal  is  a very  efficient  and  reliable  way  for 
detection  of  tool  breakage. 

4 . 5 Limitations  of  the  Cutting  Force  Signal 
for  Tool  Breakage  Detection 

Although  the  difference  per  tooth  average  method  works 
well  for  the  detection  of  tool  breakage,  measurement  of  the 
cutting  force  signal  is  difficult  in  milling.  In  the 
experimental  set-up,  an  unsatisfactory  transfer  function 
between  the  measured  force  and  the  actual  force  is 
obtained,  see  Figure  4.16.  Ideally,  the  magnitude  of 
transfer  function  at  all  frequencies  is  equal  to  one. 
Obviously,  the  dynamometer  distorts  the  cutting  force 
signal,  and  does  so  even  worse  at  high  frequencies.  It  is 
seen  that  the  dynamometer  works  well  only  under  100  Hz. 
This  is  because  the  dynamometer  with  the  workpiece  and 
with  the  angle  plate  acts  as  a system  with  several  flexible 
vibratory  modes. 

In  order  to  verify  the  distortion  of  the  cutting  force 
signal,  two  up  milling  cutting  tests  are  shown  in  Figure 
4.23.  The  cutting  parameters  are  listed  in  Table  4.5. 
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Figure  4. 21. a)  Average  force  and  difference  of  the 

average  force  in  the  1/8  immersion  of  cut 
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Figure  4.21.b)  Average  force  and  difference  of  the 

average  force  in  the  1/4  immersion  of  cut 
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Figure  4. 22. a)  Average  force  and  difference  of  the 

average  force  in  the  3/4  immersion  of  cut. 
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Figure  4.22.b)  Average  force  and  difference  of  the 

average  force  in  the  1/1  immersion  of  cut. 
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Figure  4.23.b)  Cutting  force  signals  in  the  Y and  Z 

directions  at  spindle  speed,  n = 1935  rpm. 


Table  4.5  Cutting  parameters  with  different  speeds  to 
verify  the  dynamometer  response 
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No 

spindle 

speed 

depth  of 
cut 

feed  per 
tooth 

radial 

immersion 

workpiece 

material 

rpm 

mm 

mm 

1 

600 

2.54 

0.254 

1/2 

cast  iron 

2 

1935 

It 

tl 

II 

It 

The  first  case  shows  the  cutting  test  data  in  the 
steady  state  corresponding  to  the  spindle  speed  n = 600 
rpm,  tooth  frequency  ft  = 80  Hz.  The  first  row  displays  the 
cutting  force  signal  in  the  Y direction  and  the  second  row, 
in  the  Z direction.  For  the  undamaged  cutter,  the  periodic 
variation  per  tooth  is  shown  in  Figure  4. 23. a and  Figure 
4.23.b.  For  the  damaged  cutter,  the  tool  breakage 
variation  per  revolution  is  shown  in  Figure  4.23.e  and 
Figure  4.23.f.  It  is  shown  that  the  response  of  the 
dynamometer  is  still  acceptable  at  this  tooth  frequency. 

Increasing  the  spindle  speed,  the  second  case  shows 
the  cutting  test  done  at  spindle  speed  n = 1935  rpm,  tooth 
frequency  ft  = 258  Hz.  It  becomes  very  hard  to  identify  the 
periodic  variation  in  the  Y direction  for  the  undamaged 
cutter  (see  Figure  4.23.c).  This  is  because  the  signal  of 
the  tooth  frequency  of  258  Hz  is  not  passed  by  the  low 
frequency  response  of  dynamometer.  Moreover,  the  tool 
breakage  feature  is  also  changed  (see  Figure  4.23.g). 
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Fortunately,  the  cutting  force  signal  in  the  Z direction 
(see  Figure  4.23.d  and  Figure  4.23.h)  still  truly 
illustrates  the  characteristic  changes  due  to  the  tool 
breakage  because  of  the  better  frequency  response  of  the 
dynamometer  in  the  Z direction.  Since  the  cutting  force 
signal  is  the  basis  of  the  detection  algorithm,  how  the 
distorted  cutting  force  signal  affects  the  detection 
algorithm  needs  to  be  investigated. 

Figure  4.24  shows  the  simulation  of  spectra  of  the 
resultant  cutting  force  signal  in  the  1/8,  1/4,  1/2,  and 

1/1  immersions  of  cut  during  the  steady  state  portion  of 
the  cut.  For  the  undamaged  cutter,  these  spectra  are 
similar  except  that  the  variable  component  of  the  cutting 
force  signal  is  zero  in  the  1/1  immersion  case  without  run- 
out (see  Figure  4.24.d).  For  the  damaged  cutter,  the 
components  in  the  tool  breakage  zone,  located  at  less  than 
ft  and  higher  than  the  DC  component,  clearly  stand  out. 

As  mentioned  in  Section  4.2,  sudden  and  large 
deviations  in  the  difference  of  the  average  cutting  force 
signal  are  attributed  to  the  components  in  the  tool 
breakage  zone.  For  the  damaged  cutter  in  the  1/8  and  1/4 
immersions,  the  frequency  components  inside  the  tool 
breakage  zone  spread  about  the  point  f-t/2  similar  to  a 
discrete  normal  distribution  function.  In  fact,  the  total 
of  eight  components  including  the  component  of  tooth 
frequency  stands  out  because  there  are  eight  teeth  on  the 
cutter.  With  increasing  the  radial  immersion,  the  peaks  of 
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Figure  4.24  Spectra  of  simulation  of  the  cutting  force 
in  the  1/8,  1/4,  1/2,  and  1/1  immersions 
of  cut. 
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the  tool  breakage  components  gradually  shift  close  to  fj-. 
For  the  1/1  immersion  case,  only  ft/8,  2ft/8,  and  3ft/8 
components  clearly  stand  out  in  the  tool  breakage  zone  due 
to  the  smoother  variation  of  the  tool  breakage  signal,  see 
Figure  4.4. 

In  reality,  a certain  degree  of  magnitude  change  for 
each  tool  breakage  component  is  due  to  the  cutting  force 
signal  filtered  through  the  unsatisfactory  transfer 
function  of  dynamometer.  Now  the  question  is  which 
components  are  more  important  in  their  contributions  to 
the  tool  breakage  feature,  i.e.  the  large  variation  of  the 
difference  of  average  force.  Figure  4.25.b  and  Figure 
4.25.C  show  the  peak  to  peak  ratios  of  individual  harmonic 
tool  breakage  component  (see  Fig  4. 25. a),  contributing  to 
the  average  and  the  difference  of  the  average  forces. 
Multiplying  the  ratios  as  shown  in  Figure  4.25.C  with  the 
corresponding  components  in  the  tool  breakage  zone  in  Fig 
4.24,  it  is  seen  that  the  tool  breakage  components  of  ft/8, 
2ft/8,  3ft/8,  4ft/8,  and  5ft/8  are  much  more  valuable  in 
producing  the  difference  of  the  average  force. 

To  summarize,  a reasonably  undistorted  frequency 
response  of  the  cutting  force  signal  up  to  half  of  the 
tooth  frequency  is  necessary  for  successful  detection  of 
milling  cutter  breakage  using  the  difference  per  tooth 
average  method. 
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Figure  4.25  Peak  to  peak  ratios  for  the  different 
frequency  components  contributing  to 
the  tool  breakage  feature, 
a)  force  ratios;  b)  average  force  ratios; 
c)  difference  of  the  average  force  ratios. 
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In  the  high  speed  milling,  high  tooth  frequencies  are 
common.  In  [48],  it  has  been  shown  that  the  optimal  stable 
speed  for  a high  metal  remove  rate  is  when  the  tooth 
frequency,  f^,  is  equal  to  the  natural  frequency,  f^.  The 
natural  frequency  of  the  spindle  for  the  most  significant 
mode  is  usually  in  the  range  of  500  - 600  Hz.  However,  it 
is  very  difficult  to  achieve  a reasonably  flat  response  of 
a dynamometer  with  the  experimental  set-up  at  these 
frequencies  or  even  half  of  them  [30]. 

In  addition  to  this,  the  dynamometer  also  inhibits  a 
considerable  portion  of  useful  working  space  for  milling 
operations,  and  its  cost  also  makes  it  not  suitable  for 
every  milling  machine.  Owing  to  these  disadvantages,  the 
use  of  an  alternative  signal  instead  of  the  cutting  force 
signal  should  be  explored. 


CHAPTER  5 

ANALYSIS  AND  DETECTION  OF  TOOL  BREAKAGE 
USING  VIBRATION  SIGNALS 


5 . 1 Reconstruction  of  the  Cutting  Force  Signal 
by  Using  an  Accelerometer 

In  Chapter  4,  an  efficient  detection  algorithm  for 
tool  breakage  has  been  developed.  However,  the  difficulties 
of  using  a dynamometer  to  measure  the  cutting  force  signal 
become  a major  obstacle  for  the  practical  use  of  this 
method.  This  is  especially  true  for  high  speed  milling 
where  the  tooth  frequencies  are  the  order  of  600  Hz. 
Another  way  of  obtaining  the  cutting  force  signal  must  be 
sought,  so  that  the  dynamometer  might  be  excluded. 

It  was  proposed  to  locate  two  accelerometers  in  the  X 
and  Y directions  on  the  headstock  of  milling  machine,  since 
they  are  not  expensive,  readily  available,  and  easy  to  use. 
Although  the  cutting  force  acting  on  the  cutter  will 
excite  the  vibration  on  the  headstock,  the  accelerometer 
signal  does  not  resemble  the  cutting  force  signal  at  all. 
The  task  becomes,  then,  how  to  reconstruct  the  cutting 
force  signal  from  the  accelerometer  signal. 
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The  characteristics  of  the  accelerometer  signal  can 
be  determined  by  the  transfer  function  between  the  cutting 
force  on  the  cutter  and  the  vibrational  acceleration  on 
the  headstock.  The  transfer  function,  G(o)  can  be  expressed 
as 

A(cj) 

G(o)  = (5.1) 

F(o) 

where 

A(w)  - frequency  response  of  the  excited  vibrational 
acceleration  on  the  headstock 

F(w)  - frequency  response  of  the  impulse  force  on  the 
cutter 

Considering  the  system  to  be  linear,  the  transfer 
function  may  be  mathematically  defined  as  a unique 
transformation  which  maps  the  input  of  cutting  force,  f(n) 
into  the  output  of  vibrational  acceleration,  a(n)  in  the 
discrete  time  domain.  This  is  denoted  as 

a(n)  = G[f(n)]  (5.2) 

and  is  also  expressed  as 


f(n) 


G[  ] 


a(n) 


In  other  words,  the  input  signal,  f(n)  has  been 
filtered  by  a linear  system  that  results  in  the  distorted 
signal,  a(n).  Now  the  problem  is  to  recover  the  original 
force  signal,  f(n)  by  the  transfer  function,  G[  ] and  the 
vibrational  acceleration,  a(n).  This  technique  is  generally 
referred  to  as  the  inverse  filtering  method. 
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5.1.1  Inverse  Filtering  Method 

Using  Equation  (5.1),  the  cutting  force  in  the 
frequency  domain  can  be  determined  by 
F(w)  = G"l(w)  * A(w) 

= H(w)  * A(w)  (5.3) 

where 

H((j)  - frequency  impulse  response  for  the  inverse 

filter 

Then,  the  cutting  force  signal,  f(n),  in  the  time 
domain  can  be  calculated  by  taking  the  inverse  Fourier 
transform  of  F(cij).  That  is 

f (n)  = F“^[H(w)  * A(cj)  ] 

h(n)  (*)  a(n)  (5.4) 

where 

(*)  - notation  of  convolution 

h(n)  - impulse  response  of  the  inverse  filter  in  the 
discrete  time  domain 

In  digital  signal  processing,  h(n)  has  a finite 
length  of  data,  but  a(n)  has  an  extremely  large  amount  of 
data  obtained  by  a data  acquisition  system.  A partitioning 
technique  is  required  to  convolve  this  kind  of  data  since 
the  cutting  force  signal,  f(n),  should  be  obtained  in  real- 
time for  the  purpose  of  tool  breakage  detection.  With  the 
development  of  the  Fast  Fourier  Transform  (FFT) , two 
partitioning  techniques  could  be  applied  to  solve  this 
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problem.  They  are  the  over-and-add  and  over-and-save 
methods  [49]. 

Although  this  method  could  be  implemented  in  real- 
time, a certain  time  delay  still  exists  in  waiting  for  each 
section  of  the  finite  length  of  data  input.  In  addition  to 
this,  nonlinearities  in  the  spindle  bearings  [50]  and  the 
spindle  motion  of  the  Z direction  affect  the 
characteristics  of  the  transfer  function.  Adaptivity  for 
the  identification  of  the  transfer  function  during  the 
cutting  operations  may  be  required,  and  therefore  this  is 
not  a straightforward  way  to  recover  the  cutting  force 
signal . 

5.1.2  Analog  Computer  Method 

Another  way  to  reconstruct  the  cutting  force  can  be 
realized  directly  in  the  time  domain  using  an  analog 
computer  method  instead  of  the  inverse  filtering  method.  It 
is  assumed  that  the  transfer  function  between  the  cutting 
force  and  the  vibrational  acceleration  is  due  to  a 
prominent  single  degree  of  freedom  system  with  mass  M, 
damping  C,  and  stiffness  K.  The  reconstruction  of  the 
cutting  force  signal  is  simply  to  do  the  integration  of  the 
vibrational  acceleration  and  velocity  signals.  That  is 


X = X + X * At 


(5.5) 


X = X + X * At 
M*X+C*X  + K*X  = F 


(5.6) 


(5.7) 
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where 

X - vibrational  acceleration  signal 

X - vibrational  velocity  signal 

X - vibrational  displacement  signal 

In  reality,  the  integration  must  be  made  as  a low-pass 
filter  to  avoid  the  DC  component  built  up.  This  means  that 
only  the  AC  component  of  the  cutting  force  signal  could  be 
recovered.  Furthermore,  an  inevitable  phase  distortion  is 
introduced  in  both  the  digital  and  analog  integrations. 
This  may  lead  to  considerable  errors  in  the  reconstruction 
of  the  cutting  force.  Actually  the  transfer  function  of  the 
spindle  contains  several  degrees  of  freedom,  therefore  how 
closely  the  actual  transfer  function  can  be  replaced  by  a 
prominent  single  degree  of  freedom  for  the  successful  tool 
breakage  detection  needs  to  be  further  considered. 

In  the  foregoing  discussion,  reconstruction  of  the 
cutting  force  signal  from  the  accelerometer  signal  on  the 
headstock  meets  considerable  practical  difficulties.  A more 
promising  way  using  vibration  signals  to  detect  tool 
breakage  without  the  need  to  reconstruct  the  cutting  force 
signal  has  been  studied. 

5.2  Analysis  of  Characteristic  Changes  in  Vibration 
Signals  Due  to  Tool  Breakage 

It  has  been  shown  that  the  cutting  force  signal  offers 
a very  good  feature  for  detecting  milling  cutter  breakage. 
Unfortunately,  practical  difficulties  in  obtaining  the 
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undistorted  cutting  force  signal  limit  its  applications. 
Therefore,  it  is  worthwhile  to  search  for  another  signal 
with  a tool  breakage  feature  as  good  as  the  cutting  force 
signal . 

In  this  section,  the  use  of  vibration  signals  for 
finding  tool  breakage  features  is  explored.  Two  suitable 
signals,  vibrational  acceleration  on  the  headstock  and 
vibrational  displacement  of  the  spindle,  have  been 
investigated.  Both  of  them  are  easy  to  obtain  in  milling 
machines . 

Theoretically,  it  is  assumed  here  that  the  transfer 
function  between  the  cutting  force  on  the  cutter  and  the 
vibration  on  the  headstock  or  the  vibration  of  the  spindle 
is  linear  and  has  a single  degree  of  freedom  with  mass  M = 
1 kg,  damping  ratio  Q = 0.04,  natural  frequency  fj^  = 600 
Hz,  see  Figure  5.1.  In  the  next  section,  a more  complicated 
transfer  function  with  several  degrees  of  freedom  will  be 
discussed. 

In  this  simplified  spindle  system  with  a single  degree 


of  freedom,  the  vibrational  acceleration,  x and 
displacement,  x may  be  expressed  as 

x = (F-C*x-K*x)  / M (5.8) 

X = X + X * At  (5.9) 

x=x+X*At  (5. 10) 


In  order  to  understand  the  influence  of  the  tooth 
frequency  on  the  vibration  signals  with  a broken  tooth,  a 
sequence  of  spindle  speeds,  225,  900,  2250,  4500  rpm  with 


Figure  5.1  Simulation  of  the  transfer  function  between 
the  force  and  the  vibrations. 

a)  acceleration  over  force; 

b)  displacement  over  force. 
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1/2  immersion  in  the  steady  state  milling  has  been  chosen 
for  investigation.  This  corresponds  to  a sequence  of 
ratios  of  the  natural  frequencies,  fj^/20,  fn/5,  fn/^/ 

see  Figure  5.2  and  Figure  5.3.  The  first  row  displays  the 
force  signal,  the  second  row  the  acceleration  signal,  and 
the  third  row,  the  displacement  signal. 

For  the  case  of  the  tooth  frequency  far  away  from  the 
natural  frequency,  both  the  acceleration  and  displacement 
signals  have  distinct  features  available  for  the  detection 
of  a broken  tooth  (see  Figure  5. 2. a).  It  is  seen  that  the 
force  signal  is  proportionately  reproduced  by  the 
displacement  signal.  On  the  other  hand,  the  acceleration 
signal  produces  AC  bursts  of  the  natural  frequency 
component  excited  by  the  each  tooth  and  loses  the  DC 
component.  As  a result,  increasing  the  tooth  frequency 
close  to  the  natural  frequency,  the  acceleration  signal 
dissipates  the  characteristic  changes  due  to  the  broken 
tooth.  In  contrast,  the  displacement  signal  still  retains 
the  characteristic  changes  (see  Figure  5.3.b). 

For  1/1  immersion  milling,  the  cases  of  f^  equal  to 
fn/20  and  f^  are  shown  in  Figure  5.4.  It  is  shown  that  the 
magnitude  of  the  acceleration  signal  decays  quickly  in  the 
steady  state  in  comparison  with  the  1/2  immersion  case.  In 
other  words,  the  magnitude  of  the  acceleration  signal 
attenuates  dramatically  from  the  transient  states  to  the 
steady  state  in  the  1/1  immersion  of  cut.  This  is  because 
the  acceleration  signal  cannot  pass  low  frequency 
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Figure  5.2  Simulation  of  the  force,  acceleration,  and 
displacement  signals  in  a 1/2  immersion  of 
cut  with  different  ratios  of  f-^/fn* 
a)  ft/fn  = 1/20;  b)  f^/f^  = 1/5- 
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Figure  5.3  Simulation  of  the  force,  acceleration,  and 
displacement  signals  in  a 1/2  immersion  of 
cut  with  different  ratios  of  f^/fj^. 
a)  ft/fn  = 1/2;  b)  f^/fn  = VI- 
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Figure  5.4  Simulation  of  the  force,  acceleration,  and 
displacement  signals  in  a 1/1  immersion  of 
cut  with  different  ratios  of  f^/fn* 
a)  ft/fn  = 1/20;  b)  ft/fn  = VI- 
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components  from  the  force  signal  in  the  steady  state,  see 
Figure  5.1.  This  makes  the  detection  of  tool  breakage 
using  the  acceleration  signal  more  difficult.  However,  this 
problem  does  not  exist  in  the  displacement  signal. 

In  conclusion,  the  displacement  signal  is  more  useful 
for  the  purpose  of  tool  breakage  detection.  This  is  because 
the  distinct  tool  breakage  feature  is  available  in  the 
displacement  signal  at  the  range  of  the  tooth  frequency 
below  or  equal  to  the  natural  frequency  of  the  spindle. 

5. 3 Detection  of  Tool  Breakage  Using  Displacement  Signal 

Although  the  displacement  signal  has  distinct  tool 
breakage  features,  it  is  still  difficult  to  detect  tool 
breakage  using  the  displacement  signal  directly.  This  is 
because  not  only  cut  geometry  and  run-out  change  the 
displacement,  but  also  different  tooth  frequencies 
corresponding  to  the  transfer  function  result  in  different 
magnification  ratios  for  the  displacement  signal. 
Therefore,  similar  to  the  cutting  force  signal,  it  is 
necessary  to  develop  a detection  algorithm  in  order  to 
extract  the  tool  breakage  feature  from  the  displacement 
signal. 

Again,  spectrum  analysis  is  used  to  understand  the 
tool  breakage  feature  in  the  displacement  signal.  Figure 
5.5  and  Figure  5.6  show  the  spectra  of  the  simulation  data 
in  Figure  5.2  and  Figure  5.3.  The  first  row  displays  the 
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Figure  5.5  Spectra  of  the  simulation  of  force  and 

displacement  in  the  1/2  immersion  of  cut 
with  different  ratios  of  ft/^n* 
a)  ft/fn  = V20;  b)  ft/^n  = V5. 
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Figure  5.6  Spectra  of  the  simulation  of  force  and 

displacement  in  the  1/2  immersion  of  cut 
with  different  ratios  of  f^/fn* 
a)  ft/fn  = 1/2;  b)  ft/fn  = VI- 
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force  spectrum  and  the  second  row  displays  the 

displacement  spectrum.  It  is  seen  that  the  tooth  frequency 
component  in  the  displacement  spectrum  is  equal  to  the  ^ ^ 

i % ; fo) 

magnification  ratio  of  the  transfer  function  multiplying  \ 

- 

the  corresponding  force  spectrum  component.  However,  the 

components  inside  the  tool  breakage  zone  at  frequencies 

less  than  f-j-  in  the  force  spectrum  are  reproduced  in  the 
displacement  spectrum  without  regard  to  the  changes  of 
tooth  frequency. 

Therefore,  it  is  clearly  seen  that  simply  applying  the 
difference  per  tooth  average  method  developed  in  Chapter  4 
can  extract  the  tool  breakage  feature  in  the  displacement 
signal.  Figure  5.7  and  Figure  5.8  show  simulation  results 
of  the  average  and  the  difference  of  the  average 

displacement  signals  which  correspond  to  the  displacement 
signal  in  Figure  5.2  and  Figure  5.3.  A clear  distinction 
between  the  undamaged  and  damaged  cutters  is  observed  and 
is  also  independent  of  the  tooth  frequency  changes.  Figure 
5.9  shows  the  complete  cutting  cycle  of  the  average  and  the 
difference  of  the  average  displacement  signals  for  the 
tooth  frequency  equal  to  f^-  Again,  the  feature  of  tool 
breakage  is  very  clear  not  only  in  the  steady  state  but 
also  in  the  transient  states. 

It  was  shown  by  simulation  that  the  detection  of  tool 
breakage  by  means  of  the  displacement  signal  of  the  spindle 
in  the  single  degree  of  freedom  model  is  very  successful. 

In  reality,  however,  the  spindle  system  may  contain  several 
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Figure  5.7  Simulation  of  the  displacement,  average 

displacement,  and  difference  of  the  average 
displacement  signals  with  different  ratios  of 
ft/fn-  a)  ft/fn  = V20;  b)  ft/fn  = V5. 
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Figure  5.8  Simulation  of  the  displacement,  average 

displacement,  and  difference  of  the  average 
displacement  signals  with  different  ratios  of 
ft/fn-  a)  ft/fn  = V2;  b)  ft/fn  = VI- 
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Figure  5.9  Simulation  of  the  average  displacement  and 

difference  of  the  average  displacement  signals 
in  the  complete  cutting  cycle  with  the  ratio 
of  ft/fn  = VI- 


96 


degrees  of  freedom.  Figure  5.10  shows  an  example  of  a 
transfer  function  of  the  spindle  is  shown  in  Fig  5. 10. a. 
Figure  5.10.b  shows  a "curve  fitted"  transfer  function  used 
to  obtain  the  modal  parameters  of  the  spindle  system.  The 
identified  modal  parameters  of  this  transfer  function  are 
listed  in  Table  5.1. 

Table  5.1  Modal  parameters  for  the  transfer  function 
of  the  spindle 


mode 

no 

mass 

(kg) 

damping 

ratio 

natural 
frequency  (Hz) 

1 

680.0 

0.047 

274 

2 

90.9 

0.140 

396 

3 

80.3 

0.027 

474 

4 

260.6 

0.037 

633 

5 

40.5 

0.028 

855 

In  order  to  investigate  the  detection  algorithm 
applied  to  this  spindle  system.  Figure  5.11  shows  the  two 
simulation  results  in  up  milling.  One  is  at  the  low 
spindle  speed  n = 480  rpm,  tooth  frequency  f^  = 64  Hz;  the 
other  is  at  the  high  spindle  speed,  n = 3555  rpm,  tooth 
frequency,  ft  = 474.  The  tooth  frequency  of  the  second  case 
is  exactly  equal  to  the  natural  frequency  of  the  most 
dominant  mode.  The  cutting  parameters  are  listed  in  Table 
5.2.  In  the  first  case,  the  average  displacement  signal  is 
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Figure  5.10  Spindle  system  with  multiple  degrees  of 
freedom,  a)  transfer  function  between 
the  force  and  displacement;  b)  curve 
fitted  transfer  function. 
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Figure  5. 11. a)  Simulation  of  the  average  displacement 

and  difference  of  the  average  displacement 
signals  at  spindle  speed,  n = 430  rpm. 
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Figure  S.ll.b)  Simulation  of  the  average  displacement 

and  difference  of  the  average  displacement 
signals  at  spindle  speed,  n = 3555  rpm. 
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shown  in  the  first  row.  The  corresponding  difference  of  the 
average  displacement  is  shown  in  the  second  row.  The  clear 
distinction  between  the  undamaged  and  damaged  cutters  is 
shown.  As  to  the  second  case,  the  average  and  difference  of 
the  average  displacement  signals  are  shown  in  the  third  and 
fourth  rows  respectively.  The  distinct  tool  breakage 
feature,  independent  of  the  tooth  frequency  changes,  is 
still  available  in  the  damaged  cutter  signal.  This  tells  us 
that  the  detection  algorithm  can  also  be  applied  to  the 
spindle  system  with  several  degrees  of  freedom. 

Table  5.2  Cutting  parameters  in  the  simulation  of  the 
cutting  tests  using  a spindle  system  with 
multiple  degrees  of  freedom 


No 

spindle 

depth  of 

feed  per 

radial 

workpiece 

speed 

cut 

tooth 

immersion 

material 

rpm 

mm 

mm 

1 

480 

2.54 

0.254 

1/2 

cast  iron 

2 

3555 

II 

II 

II 

II 

To  summarize,  processing  of  the  displacement  signal, 
through  simulation,  has  been  proven  to  be  feasible  for  tool 
breakage  detection.  In  addition  to  this,  it  was  shown  that 
a single  degree  of  freedom  model  is  enough  to  illustrate 
and  represent  the  tool  breakage  feature  in  the  real  spindle 
system  with  multiple  degrees  of  freedom.  Thus,  an 
experimental  verification  is  necessary. 
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5 . 4 Experimental  Set-up 

The  description  of  the  machining  center  and  the  data 
recording  system  was  given  in  Section  4.3.  For  the  face 
milling,  the  ZN3M  Carbaloy  milling  cutter  was  still  used 
with  three  kinds  of  inserts,  depending  on  the  workpiece 
material  and  axial  depth  of  cut.  The  catalog  numbers  of  the 
inserts  are  SNG634  (silicon  nitride) , SNG632  (silicon 
nitride) , and  SNE633  (carbide) . The  run-out  on  the  face 
milling  cutter  was  adjusted  within  0.0254  mm  (0.001").  For 
the  end  milling,  19.05  mm  (0.75")  diameter  carbide  end 
mill  cutters  (Robb  Jack  Co.)  with  2,  3,  and  4 flutes  of  30° 
helix  angle  was  used.  In  addition  to  this,  three  kinds  of 
workpiece  material,  GM-241  cast  iron,  aluminum  7075- 
T7351,  and  steel  4330  were  available  for  cutting  tests. 

The  displacement  signal  was  obtained  by  a multi- 
channel non-contact  measurement  device,  Accumeasure  system 
1000  (MTI  Co.).  This  system  determines  the  capacitance 
between  a capacitance  probe  with  a sensing  electrode  and  a 
ground  target.  Then,  the  capacitance  is  transmitted  to  a 
probe  amplifier  module  (AS-1023PA)  with  a -3  dB  corner 
frequency  at  5 KHz.  An  analog  signal  proportional  to  the 
probe-to-target  gap  is  generated  in  the  circuit  when  a 
resultant  variable  capacitance  appears. 

Two  capacitance  probes  were  located  in  the  X and  Y 
directions  between  the  housing  and  spindle,  close  to  its 
flanges  (see  Figure  5.12).  In  the  design  of  the  lubrication 
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system  for  the  front  bearings  of  the  spindle,  an  automatic 
mist  oil  spray  system  is  provided  in  the  machining  center. 
However,  the  lubrication  oil  flows  through  the  gap  between 
the  probe  and  the  spindle  when  the  spindle  rotates.  This 
may  cause  a slight  measurement  error  in  the  displacement 
signal . 

By  experiments,  some  high  freguency  components  were 
generated  in  the  displacement  signal,  especially  at  high 
spindle  speeds  close  to  the  natural  frequency  of  the 
spindle.  The  sampling  frequency  controlled  by  the  external 
clock  might  be  lower  than  twice  these  high  frequency 
contents  in  the  displacement  signal.  This  leads  to  the 
phenomenon  of  signal  aliasing  which  shifts  the  high 
frequency  components  into  the  low  frequency  domain,  such  as 
the  tool  breakage  zone,  after  the  signal  is  digitized. 
Therefore,  an  anti-aliasing  low-pass  filter  is  necessary 
for  high  spindle  speeds.  In  these  experiments,  a 2nd  order 
Butterworth  analog  low-pass  filter  was  designed  to  reduce 
the  aliasing  error.  Its  cutoff  frequency  is  at  1 KHz.  The 
circuit  diagram  is  shown  in  Figure  5.13.  Figure  5.14  is  a 
schematic  of  the  experimental  set-up  for  the  displacement 
measurements . 

The  spindle  run-out  was  not  taken  into  account  in 
Section  5.3  for  the  simulation  of  the  average  displacement 
signal.  In  fact,  a strong  periodic  displacement  signal  per 
revolution  is  generated  when  the  spindle  is  free  run. 
Therefore,  the  difference  per  tooth  average  method  as 
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Figure  5.12  Capacitance  probe  mounting. 


C 


Figure  5.13  Electrical  circuit  diagram  for  an 
anti-aliasing  filter. 
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Capaci'taiice  Probe 


Figure  5.14  Process  diagram  for  the  experimental  work. 
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applied  to  the  displacement  signal  needs  a little 
modification. 

Assuming  that  the  vibrational  response  of  the  spindle 
is  a linear  system,  the  actual  average  displacement  per 
tooth  can  be  calculated  by  subtracting  the  run-out  pattern 
before  cutting  from  the  measured  average  displacement  per 
tooth.  That  is 


n 

S Dx(i) 
i=l 


-"xavg 


(j)  = 


- Dxr(j) 


(5.11) 


n 


n 


S Dy(i) 

i=l 


D 


'yavg 


(j)  = 


n 


~ Dyr ( j ) 


(5.12) 


where 

Djjavg  “ actual  average  displacement  in  the  X direction 

Dyavg  “ actual  average  displacement  in  the  Y direction 

Dx  “ displacement  in  the  X direction 

Dy  - displacement  in  the  Y direction 

Djjr  “ average  run-out  displacement  pattern  in  the  X 
direction  when  the  spindle  is  free  run 
Dyj-  - average  run-out  displacement  pattern  in  the  Y 
direction  when  the  spindle  is  free  run 
j - index  for  time  increment  per  tooth  period 
n - number  of  samples  per  tooth  period 
The  average  resultant  displacement  per  tooth,  D^^g  can 
be  calculated  as 

Davg(j)  = (Dxavg^(j)  + Dyavg^ ( 3 ) ) (5. 13) 
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Then,  the  difference  of  the  average  displacement  per 
tooth,  ADavg  is  given  by 

ADavg(j)  ~ i*avg(3)  “ ^^avg(j“i)  (5.14) 

5 . 5 Experimental  Verification 

Figure  5.15  shows  the  raw  displacement  signal  in  up 
milling  with  the  spindle  run-out.  The  first  row  displays 
the  displacement  signals  in  the  X direction  and  the  second 
row,  in  the  Y direction.  The  cutting  parameters  are  listed 
in  Table  5.3.  It  is  not  easy  to  identify  which  column 
corresponds  to  the  damaged  cutter  and  which  tooth  was 
broken  due  to  the  strong  component  of  the  spindle  run-out. 
However,  applying  the  difference  per  tooth  average  method 
to  the  raw  data,  the  average  and  the  difference  of  the 
average  displacement  signals  are  shown  in  Figure  5.16.  A 
clear  distinction  is  available  between  the  undamaged  and 
damaged  cutters.  The  expanded  views  for  the  entry  and 
steady  states  are  shown  in  Figure  5.17. 

Table  5.3  Cutting  parameters  in  the  detection  of  tool 
breakage  using  displacement  signal 


spindle 

speed 

depth  of 
cut 

feed  per  radial 

tooth  immersion 

workpiece 

material 

rpm 

mm 

mm 

600 

2.54 

0.254  1/2 

cast  iron 
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Undamaged  Cutter 
* 10“5(M)  a) 

2.0 
1.0 
0.0 


Damaged  cutter 

* 10"® (M)  c) 

2.01 

J 

J 

1.0 


0.0 


Figure  5.15  Raw  displacement  signals  in  the  X and  Y 
directions. 
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Undaunaged  cutter 
* 10"® (M)  a) 


Damaged  Cutter 


* 10"® (M)  c) 


b) 


Figure  5.16  Average  displacement  and  difference  of  the 

average  displacement  signals  in  the  complete 
cutting  cycle. 
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Undamaged  Cutter 


Dzuaaged  Cutter 


* lo“^{M) 


a) 


0.0  0.4 


* 10"®(M)  e) 


t) 


(SEC) 


(SEC) 


Figure  5.17  Average  displacement  and  difference  of  the 
average  displacement  signals  in  the  entry 
and  steady  states. 
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Figure  5.18  shows  an  experimental  result  of  three 
consecutive  tool  breakages  during  milling  of  the  steel 
workpiece.  The  cutting  parameters  are  listed  in  Table  5.4. 
This  is  a very  good  example  of  how  important  it  is  to 
detect  tool  breakage  immediately. 

The  first  chipped  tooth  suddenly  happened  at  1.02  sec, 
then  the  next  one,  at  2.48  sec,  and  finally  the  third  one, 
at  3.98  sec.  If  this  damaged  cutter  was  still  in  use,  a 
catastrophic  tool  failure  would  occur  very  soon.  It  is 
shown  that  the  difference  of  the  average  displacement 
signal  clearly  indicates  where  tool  breakage  occurs.  Once 
the  first  chipped  tooth  is  detected,  the  damage  spreading 
to  the  other  teeth  could  be  prevented. 

Table  5.4  Cutting  parameters  in  up  milling  with 
three  consecutive  tool  breakages 


spindle 

speed 

depth  of 
cut 

feed  per 
tooth 

radial 

immersion 

workpiece 

material 

rpm 

mm 

mm 

600 

2.54 

0.381 

1/2 

steel 

5.6  Comparison  of  Detection  of  Tool  Breakage  Using 
Displacement  Signal  and  Force  Signal 

It  seems  that  the  average  and  the  difference  of  the 
average  displacement  signals  are  very  similar  to  those  of 
force  signals.  Figure  5.19  shows  the  average  and  the 
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* 10"5{M)  a) 


b) 


(SEC) 


Figure  5.18  Average  displacement  and  difference  of  the 
average  displacement  signals  with  three 
consecutive  tool  breakages. 


Ill 


Force 


* 10^ (N)  a) 


Displacement 

* 10“^ (M)  c) 


b)  d) 


Figure  5.19  Simulation  of  the  average  and  difference  of 
the  average  force  and  displacement  signals 
independent  of  the  tooth  frequency  changes. 


112 


difference  of  the  average  force  and  displacement  signals 
in  Figure  5.2  and  Figure  5.3.  Although  the  instantaneous 
displacement  signal  changes  depending  on  the  ratio  of 

ff/fn'  shown  that  after  averaging  and  differencing 

per  tooth,  the  displacement  signal  is  proportional  to  the 
force  signal  independent  of  the  tooth  frequency  changes. 
According  to  the  transfer  function  shown  in  Figure  5.1,  the 
ratio  between  the  average  displacement  and  force  is 

approximately  equal  to  the  DC  component  of  the  transfer 
function,  1/K,  where  K is  the  stiffness  parameter. 

K = M * (27t  * fn)2 

= 1 * (27T  * 600)2  = 1.421  * 10'7 
Actually,  the  relationship  between  the  average  force 
and  displacement  may  also  be  expressed  using  the  equations 
developed  in  Chapter  3.  For  the  face  milling  with  run-out 
free  cutter,  the  instantaneous  forces  in  the  X and  Y 
directions  may  be  expressed  as 
fx  = fT  * COS0  + fj^  * sinG 

= kg  * a^  * ft  * (sinG  * cosG  + c * sin^G)  (5.15) 

and 

fy  = -fip  * sinG  + fj^  * cosG 

= ks  * 3a  * ft  * (-sin^G  + c * sinG  * cosG)  (5.16) 

where 

kg  - specific  cutting  force 
a^  - axial  depth  of  cut 
ft  - feed  per  tooth 
c - ratio  of  f-ji/f^,  constant 
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The  corresponding  average  forces  in  the  X and  Y 
directions  are  given  by 


m 


xavg 


2n 


’ks*aa*ft* (sin0*cos9  + c*sin20)d0  (5.17) 


m 


yavg 


27T 


0, 


’ks*aa*ft* + c*sin0 *cos0 ) d0  (5.18) 


where 

xn  — number  of  the  teeth  on  the  cutter 
0g  - starting  angle  of  cut 
0g  - exit  angle  of  cut 

Then,  the  resultant  average  force  may  be  expressed  in 
the  form, 

F, 


“ (^xava^  ^yavg^)^'^^ 

F(0g,  0©/  c) 


avg  “ v^'xavg  ' '•yavg 
m*kg*aa*ft 


2tt 

m*kg*aa*ft 

27T 


F(<Psf  c) 


(5.19) 


(5.20) 


where 

F(0g/  0©/  ~ function  of  0g,  0e/  ^ 

F(0g»  “ fuJ7ction  of  0g,  c 

- swept  angle  of  cut,  i.e.  0g  ” ®s 
Finally,  the  resultant  average  displacement  may  also 

be  expressed  in  the  form, 

m*kg*aa*ft 


^avg 


27tK 


F(0s' 


(5.21) 


Figure  5.20  shows  the  experimental  comparison  between 
the  average  force  and  displacement  signals  in  a sequence  of 
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Force 

* lO^(N)  a) 

4.0  1 


Displacement 

* 10"5(M)  e) 

2.0  1 i 


0.0  3.0 

(SEC) 


Figure  5.20  Average  force  and  displacement  signals 
obtained  from  different  spindle  speeds. 
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spindle  speeds.  The  cutting  parameters  are  listed  in  Table 

5.5. 


Table  5.5  Cutting  parameters  in  up  milling  with 
force  and  displacement  sensors  active 


No 

spindle 

speed 

depth  of 
cut 

feed  per 
tooth 

radial 

immersion 

workpiece 

material 

rpm 

mm 

mm 

1 

600 

2 . 54 

0.254 

1/2 

cast  iron 

2 

1500 

II 

II 

II 

II 

3 

1845 

II 

II 

II 

II 

4 

2250 

II 

II 

II 

II 

Figure  5.21  shows  the  corresponding  difference  of  the 
average  force  and  displacement  signals.  It  is  shown  that 
the  average  or  the  difference  of  the  average  displacement 
can  be  approximately  treated  as  the  average  or  the 
difference  of  the  average  force  by  multiplying  a constant, 
K. 

In  the  next  chapter,  numerous  cutting  test  results 
will  be  presented  to  verify  this  detection  algorithm  under 
a variety  of  cutting  conditions  using  the  displacement 
signals. 
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Force  Displacement 


3.0 

0.0 

(SEC) 

(SEC) 

Figure  5.21  Difference  of  the  average  force  and 
displacement  signals  obtained  from 
different  spindle  speeds. 


CHAPTER  6 

DETECTION  OF  TOOL  BREAKAGE  USING  DISPLACEMENT  SIGNAL 
IN  A VARIETY  OF  CUTTING  CONDITIONS 


The  use  of  the  difference  per  tooth  average  method  to 
detect  tool  breakage  in  the  displacement  signal  has  been 
proven  to  be  feasible.  However,  the  applicability  of  this 
method  needs  to  be  verified  to  various  cutting  conditions. 
A lot  of  experimental  results  using  different  cutting 
parameters  are  shown  in  this  chapter. 

6 . 1 Spindle  Speed 

Figure  6.1  shows  the  transfer  function  between  the 
force  on  the  cutter  and  the  displacement  on  the  spindle  in 
the  Y and  X directions  in  the  experimental  set-up  as  shown 
in  Figure  5.12.  There  are  several  modes  with  resonant  peaks 
in  both  directions.  Within  the  range  of  interest,  it 
contains  two  resonant  peaks;  at  262  Hz,  476  Hz  in  the  Y 
direction  and  at  244  Hz,  450  Hz  in  the  X direction.  By 
experiments,  the  transfer  function  is  changed  little  by  the 
spindle  motion  of  the  Z axis. 
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* 10“®  (M/N)  a) 


Figure  6.1  Transfer  function  between  the  force  on  the 
cutter  and  displacement  of  the  spindle, 
a)  Y direction;  b)  X direction. 
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Different  spindle  speeds  have  their  tooth  frequencies 
corresponding  to  the  magnification  ratios  in  Figure  6.1. 
In  up  milling  of  cast  iron  using  the  face  milling  cutter, 
the  displacement  signal  in  the  Y direction  is  the  major 
component  in  comparison  to  that  in  the  X direction. 
Therefore,  a sequence  of  cutting  tests  was  done,  and  the 
locations  of  tooth  frequencies  corresponded  to  the  marks  in 
the  transfer  function  of  the  Y direction.  Figure  6.2  shows 
the  average  displacement  signal,  and  the  difference  of  the 
average  displacement  signal  is  shown  in  Figure  6.3.  The 
cutting  parameters  are  listed  in  Table  6.1. 

Table  6.1  Cutting  parameters  in  up  milling  with 
different  spindle  speeds 


No 

spindle 

depth  of 

feed  per 

tooth 

radial 

speed 

cut 

tooth 

frequency 

immersion 

rpm 

mm 

mm 

Hz 

1 

600 

1.27 

0.254 

80 

1/2 

2 

1425 

II 

II 

190 

II 

3 

1965 

II 

II 

262 

II 

4 

2160 

II 

II 

290 

II 

5 

3000 

II 

It 

400 

II 

6 

3570 

II 

If 

476 

II 

In  down  milling,  the  displacement  signal  in  the  X 
direction  is  the  major  component.  Therefore,  another 
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C) 

i) 
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• 

d) 

j) 

0. 

10. 

5) 

5. 

0.0  3.333 

0.0  3.333 

e) 

k) 
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10. 

0.0  2.40 

0.0  2.40 

6) 

5. 

f) 

1) 
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°*  0.0  2.017  0.0  2.017 
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Figure  6.2  Average  displacement  signal  in  up  milling  with 
different  spindle  speeds. 
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1) 
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Figure  6.3  Difference  of  the  average  displacement  signal 
in  up  milling  with  different  spindle  speeds. 
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sequence  of  cutting  tests  was  done,  see  Figure  6.4  and 
Figure  6.5.  The  tooth  frequencies  correspond  to  the  marks 
in  the  transfer  function  of  the  X direction.  The  cutting 
parameters  are  listed  in  Table  6.2. 

Table  6.2  Cutting  parameters  in  down  milling  with 
different  spindle  speeds. 


No 

spindle 

depth  of 

feed  per 

tooth 

radial 

speed 

cut 

tooth 

frequency 

immersion 

rpm 

mm 

mm 

Hz 

1 

600 

1.27 

0.254 

80 

1/2 

2 

1830 

II 

II 

244 

II 

3 

1920 

II 

II 

256 

II 

4 

3000 

II 

II 

400 

It 

5 

3375 

tf 

II 

450 

II 

It  is  seen  that  the  displacement  signal  can  detect 
tool  breakage  quite  well,  even  at  a tooth  frequency  equal 
to  the  most  dominant  resonant  mode.  As  to  the  most  dominant 
mode  cases  in  up  and  down  millings,  the  difference  signal 
becomes  a little  bigger  for  the  undamaged  cutter  due  to  the 
large  vibration  and  high  frequency  noise  (see  Figure  6.2.f 
and  Figure  6.4.e)  Therefore,  it  may  make  the  detection 
margin  smaller  between  undamaged  and  damage  cutters. 
However,  experiments  demonstrated  that  this  could  be 
further  improved  by  using  a high  order  anti-aliasing  filter 
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Figure  6.4  Average  displacement  signal  in  down  milling 
with  different  spindle  speeds. 
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Figure  6.5  Difference  of  the  average  displacement  signal 
in  down  milling  with  different  spindle  speeds. 
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so  that  the  margin  increases. 

As  a result,  the  tool  breakage  feature  is  almost 
independent  of  the  change  of  spindle  speeds  as  shown  in  the 
simulation  results  in  Chapter  5.  This  is  because  the 
freguency  response  of  the  transfer  function  is  reasonably 
flat  up  to  half  or  over  half  of  the  dominant  natural 
frequency  (see  the  explanations  in  Sections  4.5  and  5.6) . 

6 . 2 Workpiece  Material 

In  Equation  (5.21),  the  average  displacement  signal 
is  proportional  to  the  cutting  stiffness.  Therefore,  the 
different  workpiece  materials  will  have  different  levels  of 
detection  signal  in  the  difference  of  the  average 
displacement,  even  using  the  same  other  cutting  parameters. 

Figure  6.6  shows  the  average  displacement  signal 
obtained  for  three  kinds  of  workpiece  material  in  up 
milling.  The  cutting  parameters  are  listed  in  Table  6.3. 
The  first  row  displays  the  signal  obtained  from  the  steel 
workpiece.  The  second  row  displays  the  signal  obtained  from 
the  cast  iron  workpiece  and  the  third  row,  the  aluminum 
workpiece.  It  is  clearly  shown  that  the  order  of  cutting 
stiffness  is  steel,  then  cast  iron,  and  then  aluminum. 

The  corresponding  difference  of  the  average 
displacement  signal  is  shown  in  Figure  6.7.  This  tells  us 
that  the  tool  breakage  signal  in  the  difference  of  the 
average  displacement  signal  is  proportional  to  the  cutting 
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(SBC) 


(SEC) 


Figure  6.6  Average  displacement  signal  obtained  from  the 
three  kinds  of  workpiece  material;  steel,  cast 
iron,  aluminum. 
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Figure  6.7  Difference  of  the  average  displacement  signal 
obtained  from  the  three  kinds  of  workpiece 
material:  steel,  cast  iron,  aluminum. 
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stiffness  of  the  workpiece  material. 


Table  6.3  Cutting  parameters  in  the  cutting  tests  using 
different  workpiece  materials 


No 

spindle 

speed 

depth  of 
cut 

feed  per 
tooth 

radial 

immersion 

workpiece 

material 

rpm 

mm 

mm 

1 

600 

1.27 

0.254 

1/2 

steel 

2 

fl 

II 

II 

II 

cast  iron 

3 

II 

II 

II 

II 

aluminum 

6 . 3 Axial  Depth  of  Cut 

The  changes  of  axial  depth  of  cut  (ADC)  can  be  divided 
by  two  kinds:  one  is  sudden  change  of  ADC,  and  the  other 
is  gradual  change  of  ADC.  Figure  6.8  shows  the  geometries 
of  sudden  and  gradual  changes  of  ADC  in  the  cutting  tests. 
The  cutting  parameters  are  listed  in  Table  6.4.  Figure  6.9 
shows  the  average  displacement  signal,  and  Figure  6.10 
shows  the  difference  of  the  average  displacement  signal. 

Except  in  the  transient  states,  the  changes  of  the 
average  and  the  difference  of  the  average  signals  are 
actually  followed  the  changes  in  the  axial  depths  of  cut. 
This  tells  us  that  the  tool  breakage  signal  in  the 
difference  of  the  average  displacement  is  proportional  to 
the  axial  depth  of  cut. 
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Figure  6.8  Cut  geometry  with  the  varying  axial  depths 
of  cut. 
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Figure  6.9  Average  displacement  signal  with  the  changes 
of  the  axial  depths  of  cut. 
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Figure  6.10  Difference  of  the  average  displacement  signal 
with  the  changes  of  the  axial  depths  of  cut. 
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Table  6.4  Cutting  parameters  in  the  cutting  tests  with 
different  axial  depths  of  cut 


No 

spindle 

speed 

depth  of 
cut 

feed  per 
tooth 

radial 

immersion 

workpiece 

material 

rpm 

mm 

mm 

1 

600 

* 

0.254 

1/2 

cast  iron 

2 

It 

* 

II 

II 

II 

3 

II 

* 

II 

II 

II 

* see  Figure  6.8 


6 . 4 Feedrate 


Feedrate  is  also  a very  important  cutting  parameter 
relating  to  tool  life.  Figure  6.11  shows  the  feedrate 
changes  in  the  cutting  tests.  Figure  6.12  and  Figure  6.13 
show  the  average  and  the  difference  of  the  average  signals 
in  the  varying  feedrates.  The  cutting  parameters  are  listed 
in  Table  6.5. 

Table  6 . 5 Cutting  parameters  in  the  cutting  tests  with 
different  feedrates 


No 

spindle 

speed 

depth  of 
cut 

feed  per 
tooth 

radial 

immersion 

workpiece 

material 

rpm 

mm 

mm 

1 

600 

1.27 

* 

1/2 

cast  iron 

2 

It 

tt 

* 

It 

tt 

3 

It 

tt 

* 

It 

tt 

* see  Figure  6.11 
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Figure  6.11  Cut  geometry  with  the  varying  feedrates. 
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Figure  6.12  Average  displacement  signal  with  varying 
feedrates . 
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Figure  6.13  Difference  of  the  average  displacement  signal 
with  varying  feedrates. 
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It  is  shown  that  the  changes  of  feedrate  affect  the 
average  and  the  difference  of  the  average  displacement 
signals  proportionally. 


6 . 5 Radial  Immersion  of  Cut 


In  milling  operations,  the  radial  immersion  of  cut  is 
changed  only  in  the  transient  state  because  of  the 
variation  of  the  swept  angle  of  cut.  Figure  6.14  shows  the 
average  displacement  signal  in  the  1/8,  1/4,  1/2,  3/4,  1/1 
radial  immersions  up  milling.  The  corresponding  difference 
of  the  average  signal  is  shown  in  Figure  6.15.  The  cutting 
parameters  are  listed  in  Table  6.6. 

Table  6.6  Cutting  parameters  in  the  cutting  tests  with 
different  radial  immersions  of  cut 


No 

spindle 

speed 

depth  of 
cut 

feed  per 
tooth 

radial 

immersion 

workpiece 

material 

rpm 

mm 

mm 

1 

600 

1.27 

0.254 

1/8 

steel 

2 

II 

II 

II 

1/4 

II 

3 

II 

II 

II 

1/2 

II 

4 

II 

II 

II 

3/4 

It 

5 

II 

II 

It 

1/1 

It 

6 

It 

It 

It 

* 

cast  iron 

7 

II 

II 

II 

* 

II 

* see  Figure  6.16 
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Figure  6.14  Average  displacement  signal  in  the  1/8,  1/4, 
1/2,  3/4,  and  1/1  immersions  of  cut. 
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Figure  6.15  Difference  of  the  average  displacement  signal 
in  the  1/8,  1/4,  1/2,  3/4,  and  1/1  immersions 
of  cut. 
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It  is  shown  that  the  DC  component  is  almost 
proportional  to  the  radial  immersion  of  cut,  but  the  AC 
component  of  the  tool  breakage  signal  in  the  difference  of 
the  average  displacement  is  not.  The  tool  breakage  signal 
reaches  the  maximum  value  in  the  1/2  immersion  of  cut. 
Then,  it  decreases  as  the  radial  immersion  increases. 

In  addition  to  this,  two  special  cut  geometries  with 
gradual  change  of  radial  immersion  of  cut  are  shown  in 
Figure  6.16.  Figure  6.17  shows  the  cutting  test  results. 
The  cutting  parameters  are  also  listed  in  Table  6.6.  It  is 
shown  that  clear  distinction  of  the  tool  breakage  feature 
is  available  even  in  varying  immersions  of  cut. 

6.6  Interrupted  Cuttings 

In  general,  the  boundary  of  the  workpiece  is  a 
rectangle  whose  sides  are  parallel  or  perpendicular  to  the 
feed  direction.  As  the  cutter  enters  into  the  workpiece, 
there  are  three  common  cutting  states:  entry,  steady,  and 
exit.  If  the  workpiece  contains  holes,  bores,  slots, 
pockets,  or  curvilinear  boundaries  interrupted  cuttings 
happen,  which  change  start  and  exit  angles  of  cut 
dependent  on  the  geometry  of  workpiece.  As  a result,  the 
shape  of  the  displacement  signal  may  be  different. 
Therefore,  it  is  necessary  to  investigate  the  influence  of 
interrupted  cuttings  in  the  average  and  the  difference  of 
the  average  displacement  signals. 
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Figure  6.16  Cut  Geometry  for  the  gradual  change  of 
the  radial  immersion  of  cut. 


141 


Damaged  Cutter 


e) 


5. 


d) 


0 — mmm 


J- 


0.0 


6.0 

(SEC) 


12.0  0.0 


12.0 


Figure  6.17  Average  displacement  and  difference  of 

the  average  displacement  signals  with  the 
gradual  change  of  radial  immersion  of  cut. 
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Three  special  cases  of  interrupted  cutting  are 
represented  in  this  section.  The  face  milling  cutter  was 
used  to  mill  over  a slot,  a hole,  circular  exit  and  entry 
boundaries.  The  cutting  parameters  are  listed  in  Table  6.7. 

Table  6.7  Cutting  parameters  in  the  interrupted  cuttings 


No 

spindle 

speed 

depth  of 
cut 

feed  per 
tooth 

radial 

immersion 

workpiece 

material 

rpm 

mm 

mm 

1 

600 

1.27 

0.254 

3/4 

cast  iron 

2 

II 

II 

II 

1/2 

II 

3 

II 

2.54 

II 

1/2 

II 

6.6.1  Milling  over  a Slot 

The  geometry  of  the  workpiece  is  shown  in  Figure  6.18. 
The  cutting  test  results  are  shown  in  Figure  6.19.  The 
first  row  displays  the  average  displacement  signal  and  the 
second  row,  the  difference  of  average  displacement  signal. 
It  is  shown  that  the  difference  per  tooth  average  method  is 
able  to  remove  the  transient  effect  of  milling  over  a slot. 

6.6.2  Milling  over  a Hole 

Figure  6.20  shows  the  geometry  of  the  workpiece  with  a 
hole.  The  cutting  test  results  are  shown  in  Figure  6.21. 
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Figure  6.18  Cut  geometry  of  the  slot. 
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Figure  6.19  Average  displacement  and  difference  of 
the  average  displacement  signals  during 
milling  over  a slot. 
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Figure  6.20  Cut  geometry  of  the  hole. 
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Figure  6.21  Average  displacement  and  difference  of 
the  average  displacement  signals  during 
milling  over  a hole. 
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Again,  the  first  row  displays  the  average  displacement 
signal  and  the  second  row,  the  difference  of  the  average 
displacement  signal.  In  this  case,  the  average  displacement 
signal  decreases  when  the  cutter  starts  milling  over  the 
hole.  Then,  the  signal  gradually  increases  until  the  cutter 
has  passed  the  hole.  Nevertheless,  it  is  efficient  to 
remove  the  influence  of  this  transient  and  keep  the  tool 
breakage  feature  by  using  the  difference  per  tooth  average 
method  (see  Figure  6.21.b  and  Figure  6.21.d). 

6.6.3  Milling  over  Circular  Exit  and  Entry  Boundaries 

Figure  6.22  shows  circular  exit  and  entry  boundaries. 
The  cutting  test  results  are  shown  in  Figure  6.23.  Using 
the  undamaged  cutter,  the  average  displacement  signal  is 
displayed  in  the  first  row  and  the  difference  of  the 
average  displacement  is  display  in  the  second  row. 

Due  to  the  sudden  change  of  the  workpiece  boundary,  a 
peak  of  the  difference  signal  is  generated  even  for  the 
undamaged  cutter,  which  may  be  higher  than  the  tool 
breakage  signal  produced  by  the  damaged  cutter.  However, 
this  is  only  a one  shot  event.  A special  detection  strategy 
(Section  7.1)  is  necessary  to  avoid  a false  alarm  caused  by 
these  transients. 
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Figure  6.22  Cut  geometry  with  the  circular  exit  and  entry 
boundaries . 
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Figure  6.23  Average  displacement  and  difference  of 
the  average  displacement  signals  during 
milling  over  circular  boundaries. 
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6.7  Corner  Slotting 

In  milling  operations,  there  are  two  common  corner 
slottings,  square  corner  slotting  and  corner  slotting  with 
a radius.  The  transient  response  is  much  stronger  in  the 
square  corner  slotting  because  the  cutter  switches  from 
full  immersion  to  half  immersion  at  the  corner,  then  the 
cutter  reaches  full  immersion  again,  see  Figure  6.24. 

A three  fluted  undamaged  end  mill  cutter  (diameter  = 
19.05  mm)  was  used  in  the  cutting  test.  The  cutting 
parameters  are  listed  in  Table  6.8. 


Table  6 . 8 Cutting  parameters  in  the  corner  slotting 


spindle 

speed 

depth  of 
cut 

feed  per 
tooth 

radial 

immersion 

workpiece 

material 

rpm 

mm 

mm 

600 

2.54 

0.254 

1/1 

cast  iron 

- 

Figure  6. 25. a shows  the  average  signal  and  Figure 
6.25.b,  the  difference  of  the  average  signal.  Similar  to 
the  milling  over  circular  boundaries,  either  a negative 
peak  or  a positive  peak  with  a small  time  delay  is 
generated  during  the  sudden  changes  of  immersions  at  the 
corner.  However,  this  is  still  a one  shot  event,  unlike  the 
tool  breakage  event  which  repeats  once  per  revolution  after 
a broken  tooth  exists  and  has  a minus-plus  or  plus-minus 


sequence  of  strong  signal. 
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Figure 


Figure  6. 


6.24  Cut  geometry  of  the  corner  slotting. 
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25  Average  displacement  and  difference  of 
the  average  displacement  signals  during 
slotting  the  square  corner. 
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6.8  Different  Flutes  of  Cutters 

As  shown  in  Equation  (5.21),  the  average  displacement 
signal  is  also  proportional  to  the  number  of  teeth  on  the 
cutter.  Figure  6.26  shows  slotting  results  using  the 
undamaged  end  mill  cutters  with  different  flutes.  The  first 
row  displays  the  signal  obtained  from  a 4 fluted  cutter. 
The  second  row  displays  the  signal  obtained  from  a 3 fluted 
cutter  and  the  third  row,  a 2 fluted  cutter.  The  cutting 
parameters  are  listed  in  Table  6.9. 

Table  6.9  Cutting  parameters  in  the  cutting  tests  with 
different  flutes  of  cutters 


No 

spindle 

speed 

depth  of 
cut 

feed  per 
tooth 

number  of 
teeth 

workpiece 

material 

rpm 

mm 

mm 

1 

400 

1.27 

0.254 

4 

cast  iron 

2 

600 

II 

II 

3 

II 

3 

800 

II 

It 

2 

II 

The  figures  of  the  average  signal  are  shown  in  the 
left-hand  side,  and  the  figures  of  the  difference  of  ;the 
average  signal  are  shown  in  the  right-hand  side.  It  is 
shown  that  the  DC  component  of  the  average  signal  is 
proportional  to  the  number  of  flutes  on  the  cutter.  The 
difference  of  the  average  signal  varies  depending  not  only 
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Figure  6.26  Average  displacement  and  difference  of  the 
average  displacement  signals  in  the  end 
milling  with  4,  3,  and  2 flutes  of  cutters. 
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on  the  magnitude  of  run-out  on  the  end  mill  cutter  but 
also  on  the  number  of  teeth,  which  is  not  a linear 
relationship. 


6 . 9 Percentage  of  Tool  Breakage 

Chipping  occurs  when  the  cutting  edge  is  broken  off 
only  partly  rather  than  being  totally  broken.  This  makes 
the  cutting  edge  ragged  and  inefficient.  As  a result,  the 
total  breaking  of  cutting  edge  usually  follows  very 
quickly.  Therefore,  detection  of  chipping  is  also 
important. 

Figure  6.27  shows  the  up  milling  cutting  test  results 
with  different  percentages  of  tool  breakage.  The  cutting 
parameters  are  shown  in  Table  6.10. 

Table  6.10  Cutting  parameters  in  the  cutting  tests  with 
different  percentages  of  broken  edge 


spindle 

speed 

depth  of 
cut 

feed  per 
tooth 

radial 

immersion 

workpiece 

material 

rpm 

mm 

mm 

600 

1.27 

0.254 

1/2 

cast  iron 

A face  milling  cutter  with  eight  unbroken  inserts  was 
used,  see  the  first  row.  Then,  one  of  the  inserts  was 
replaced  with  a chipped  insert  with  about  60%  of  the 
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Figure  6.27  Average  displacement  and  difference  of  the^ 
average  displacement  signals  with  different 
percentages  of  a broken  cutting  edge. 
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cutting  edge  broken  off,  see  the  second  row.  Finally,  a 
totally  broken  insert  was  used,  see  the  third  row.  Again, 
the  figures  of  the  left-hand  side  show  the  average  signal; 
the  figures  of  the  right-hand  side,  the  difference  of  the 
average  signal. 

It  seems  that  the  tool  breakage  signal  in  the 
difference  of  the  average  displacement  is  approximately 
proportional  to  the  percentage  of  broken  cutting  edge. 
Therefore,  the  difference  of  the  average  displacement 
signal  with  a proper  threshold  value  can  also  be  used  to 
detect  chipping  in  milling. 


CHAPTER  7 

STRATEGIES  FOR  DETECTION  OF  TOOL  BREAKAGE  UNDER 
VARYING  CUTTING  CONDITIONS 


In  Chapter  6,  it  was  clearly  shown  that  the  difference 
of  the  average  displacement  signal  is  effective  in  the 
recognition  of  tool  breakage.  However,  a knowledge  of 
allowable  level  of  the  difference  of  the  average 
displacement  is  necessary  to  distinguish  between  the 
undamaged  and  damaged  cutters.  Furthermore,  an  algorithm  is 
required  to  separate  tool  breakage  from  the  sudden 
transient  changes  like  corner  slotting  or  milling  over 
circular  exit  and  entry  boundaries. 


7.1  Pattern  Analysis 

In  Chapter  6,  it  was  shown  that  the  difference  of  the 
average  displacement  can  easily  remove  the  effect  of  slow 
varying  transients.  These  slow  varying  transients  include 
milling  during  entry,  exit,  slot,  hole,  or  circular 
interpolation,  etc.  However,  a problem  occurs  in  the  sudden 
transient  changes.  A sudden  one  shot  of  difference  signal 
is  generated.  In  order  to  avoid  the  occurrence  of  a false 
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alarm,  a special  algorithm  to  separate  tool  breakage  from 
all  kinds  of  transients  has  been  developed  [29].  This  is 
possible  because  the  pattern  of  sudden  transient  changes 
is  different  from  that  of  tool  breakage  or  chipping. 

It  is  shown  throughout  that  the  difference  of  the 
average  displacement  signal  decreases  then  immediately 
increases,  or  increases  then  immediately  decreases,  at  the 
presence  of  tool  breakage.  However,  the  difference  of  the 
average  displacement  signal  in  transients  without  tool 
breakage  is  either  positive  or  negative.  That  is,  a one 
shot  event. 

In  order  to  distinguish  them,  two  threshold  values  are 
required.  One  is  on  the  positive  side;  the  other  is  on  the 
negative  side.  If  the  difference  signal  is  greater  than  the 
positive  threshold  value,  the  flag  is  equal  to  1.  When  the 
difference  signal  is  smaller  than  the  negative  threshold 
value,  the  flag  is  equal  to  -1.  When  the  difference  signal 
is  within  the  band  of  the  two  threshold  values,  the  flag  is 
equal  to  0.  Tool  breakage  only  occurs  at  the  changes  of  the 
flag  from  -1  to  1 or  from  1 to  -1. 

7.2  Threshold  Setting 

For  a practical  sensor  system,  a high  sensitivity  to 
tool  breakage  and  a low  sensitivity  to  the  variation  of 
cutting  parameters  such  as  speed,  axial  depth  of  cut, 
radial  immersion  of  cut,  feed,  and  workpiece  material  are 
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required.  Therefore,  how  to  set  an  alarm  threshold  on  the 
positive  and  negative  sides  needs  to  be  investigated. 

In  the  processing  of  the  displacement  signal,  the 
average  displacement  per  tooth  is  calculated  and 
differenced  to  remove  the  regular  periodic  and 
deterministic  trends.  Generally  speaking,  small  residuals 
of  the  difference  of  the  average  displacement  are  left  for 
the  undamaged  cutter.  Hence,  a minimum  allowable  level  of 
this  difference  signal,  which  is  called  the  threshold 
value,  has  to  be  determined. 

This  small  residual  signal  is  basically  composed  of 
two  components.  One  is  the  component  of  process  noise;  the 
other  is  the  component  of  cutter  run-out.  The  amplitude  of 
the  noise  signal  can  only  be  determined  by  experiments.  In 
practice,  the  signal-to-noise  ratio  is  large  and  therefore 
the  residual  signal  can  be  mostly  attributed  to  the  cutter 
run-out.  Run-out  on  the  cutter,  especially  radial  throw, 
also  produces  the  same  pattern  in  the  difference  signal  as 
a broken  tooth.  If  the  radial  throw  becomes  severe,  it  is 
as  equally  harmful  as  tool  breakage  because  of  the  overload 
on  the  following  tooth.  It  is  then  conceivable  to  set  an 
alarm  threshold  equal  to  the  summation  of  the  noise  signal 
and  the  difference  signal  produced  by  the  tolerable  amount 
of  run-out  on  the  cutter. 

Using  the  run-out  model  developed  in  Chapter  3,  the 
difference  signal  produced  by  the  cutter  run-out  can  be 
calculated.  However,  it  is  necessary  to  know  the  cutting 
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parameters  first.  In  reality,  the  axial  depth  of  cut  and 
radial  immersion  of  cut  may  vary  in  the  cutting  process. 
Therefore,  the  development  of  strategies  for  the  real-time 
identification  of  the  cutting  geometry  is  required  so  that 
an  analytical  estimate  of  the  expected  threshold  may  be 
made. 

Unfortunately,  up  to  now,  an  efficient  method  for  the 
identification  of  the  axial  depth  of  cut  or  radial 
immersion  of  cut  is  not  available.  A similar 
identification  strategy  using  the  average  cutting  force  has 
been  presented  [51].  It  has  been  shown  that  it  is  only 
possible  to  identify  the  axial  and  radial  depths  of  cut  in 
the  steady  state.  Even  though  this  strategy  can  be  applied 
to  the  displacement  signal,  it  is  still  not  satisfactory 
due  to  the  poor  identification  of  transients. 

An  alternative  way,  which  is  a normalization  of  the 
difference  of  the  average  displacement  signal,  has  been 
proposed  [4].  It  is  proposed  to  obtain  a ratio  of  the 
difference  signal  to  a reference  signal  so  that  the  ratio 
may  be  independent  of  cutting  parameters.  In  addition,  the 
tool  breakage  feature  in  the  spectrum  is  also  available 
for  the  calculation  of  the  ratio,  see  Figure  4.24. 
Considerable  efforts  have  been  devoted  to  this  approach.  So 
far,  it  is  still  not  fully  developed. 

The  reference  signals,  either  DC  or  AC,  such  as  the 
moving  average  displacement  signal  or  the  tooth  frequency 
signal  in  the  instantaneous  displacement  signal,  have  been 
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considered  either  in  the  time  or  in  the  frequency  domains. 
For  example,  the  moving  average  displacement  signal  may  be 
expressed  as 

1 +0) 

^mavg(j)  = ^ ^avg(^)  (7-1) 

2oj  + l i=-w 

where 

Dmavg  “ moving  average  displacement 
2(j+l  - length  of  the  window  for  moving  average 
j - index  for  time  increment  per  tooth  period 

The  ratio  eliminates  the  influence  of  the  axial  depth 
of  cut,  feedrate,  workpiece  material,  and  some  cutter 
constants  because  the  reference  signal  is  also  proportional 
to  these  cutting  parameters.  However,  they  all  suffer  the 
same  problem,  which  is  that  the  ratio  is  extremely 
dependent  on  the  radial  immersion  of  cut. 

Figure  7.1  shows  a simulation  of  the  ratio  of  the 
peak  difference  signal  to  the  moving  average  displacement 
signal  per  revolution  versus  the  swept  angle  using  a face 
milling  cutter  with  eight  teeth.  The  line  of  tool  breakage 
represents  a cutting  edge  totally  broken;  the  line  of  run- 
out represents  the  radial  throw  amounting  to  one  tenth  of 
the  feed  per  tooth.  Using  Equation  (3-39),  the  two  lines 
intersects  at  = 0.  That  is 

r^.Q^  ~ ^i  “ O'lft  “ ftsinG 

and 

0 = sin”^0.1  = 5.74* 
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That  tells  us  that  radial  throw  and  tool  breakage  are 
indistinguishable  when  the  swept  angle  is  smaller  than 
5.74°.  As  shown  in  Figure  7.1,  the  ratio  of  run-out 
increases  in  the  smaller  radial  immersion  of  cut  for  the 
undamaged  cutter.  However,  the  ratio  of  tool  breakage 
decreases  quickly  in  the  larger  radial  immersion  of  cut, 
even  for  the  damaged  cutter.  In  Section  6.5,  the 
experimental  results  also  show  the  same  conclusion.  In 
other  words,  an  on-line  identification  of  the  extent  of  the 
radial  immersion  of  cut  is  still  necessary  so  that  the 
alarm  threshold  can  be  adaptively  varied  as  the  amount  of 
immersion  changes,  see  Figure  7.2.  The  solution  of  this 
approach  still  depends  on  the  development  of  an 
identification  strategy.  Furthermore,  the  signal-to-noise 
ratio  is  decreased  in  small  radial  immersions  of  cut 
because  the  moving  average  displacement  is  relatively 
small.  Therefore,  the  possibility  of  false  alarm  increases 
even,  with  the  identification  strategy,  in  small  radial 
immersions  of  cut.  This  makes  the  scheme  difficult  to 
implement. 

Figure  7.2  shows  an  example  of  the  ratio  calculated 
for  the  complete  cycle  of  cutting  tests.  The  cutting 
parameters  are  same  as  listed  in  Table  6.10.  The  first  row 
displays  the  moving  average  displacement  signal  which  is 
independent  of  tool  breakage.  The  second  row  displays  the 
difference  signals  per  tooth;  the  third  row,  the  ratio.  It 
is  seen  that  the  ratio  increases  in  the  entry  and  more  so 
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Figure  7.2  Moving  average  and  difference  of  the  average 
displacement  signals  with  the  normalization 
ratio  of  the  difference  of  the  average 
displacement  signal. 
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in  the  exit,  even  for  the  undamaged  cutter.  This  is 
because  the  entry  swept  angle  of  cut  changes  from  0°  to  90° 
and  the  exit  swept  angle  of  cut  changes  from  90°  to  0 . 

Another  promising  way  has  been  considered  [52].  It 
utilizes  a geometric  modelling  system  to  predict  the 
cutting  torque  in  face  and  end  milling  operations.  In  this 
system,  a geometric  simulation  of  the  cut  is  made  and  then 
the  cutting  parameters  may  be  estimated  from  the  input  data 
and  the  shape  of  the  removed  segment.  Therefore,  using  the 
estimated  cutting  parameters,  the  proper  threshold  values 
could  be  set  under  varying  cutting  conditions.  This  should 

be  dealt  with  in  future  research. 

In  recent  years,  commercial  tool  condition  monitoring 
systems  use  a "teaching”  method  to  obtain  the  threshold 
under  varying  cutting  conditions,  for  instance  Promess 
sensors  [53],  TCM  sensors  [54],  and  Feed  Force  Sensor  [55]. 
Even  though  this  teaching  method  has  some  disadvantages 
[56],  these  system  are  still  widely  accepted  by  the 

industry. 

This  teaching  method  can  also  be  applied  to  the 
difference  of  the  average  displacement  signal  very  well 
because  the  difference  signal  provides  the  distinct  tool 
breakage  feature.  The  teaching  method  applied  to  the 
displacement  sensor  can  be  explained  as  follows:  at  the 
beginning  of  the  milling  process,  the  system  is  set  to  a 
"teach  in"  mode,  and  a protection  device  is  turned  on.  The 
monitoring  system  measures  the  whole  cycle  of  the 
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difference  signal  and  then  calculates  limit  patterns  for 
tool  breakage  at  different  time  intervals.  In  the 

following  cuts,  the  monitoring  system  is  set  to  a "pattern 
matching"  mode.  The  limit  patterns  which  have  controllable 
factors  are  compared  with  the  new  measurement  of  the 

difference  signal.  If  it  is  over  a certain  level  and  the 
change  of  flag  is  also  from  -1  to  1 or  1 to  -1,  tool 
breakage  is  assumed. 

This  teaching  method  is  free  of  modelling  and  easy  to 
implement.  However,  a trial  cut  is  required  for 
initialization.  Using  this  teaching  methodology,  a tool 

breakage  monitoring  system  has  been  implemented  at  the 

Machine  Tool  Laboratory,  University  of  Florida.  It 
demonstrates  a high  sensitivity  and  reliability  for  the 
detection  of  tool  breakage  in  milling. 


CHAPTER  8 

CONCLUSIONS  AND  FUTURE  WORK 


A tool  diagnostic  system  has  been  developed  to  sense 
cutter  breakage  by  monitoring  the  pattern  changes  of  the 
cutting  force.  The  variation  per  tooth  is  filtered  out  by 
averaging  the  cutting  force  per  tooth.  The  pattern  of  tool 
breakage  is  extracted  by  differencing  the  average  cutting 
force. 

Because  of  the  practical  difficulties  for  measurement 
of  the  cutting  force  in  milling,  the  use  of  vibration 
signals  was  explored.  It  was  shown  the  displacement  signal 
of  the  spindle  maintains  the  same  tool  breakage  pattern  as 
the  cutting  force  signal  under  all  kinds  of  cutting 
conditions.  Finally,  the  setting  of  threshold  values  was 
also  investigated  to  make  the  system  independent  of  changes 
in  cutting  conditions. 

Further  research  should  be  focused  on  the  development 
of  a real-time  strategy  for  the  identification  of  cut 
geometry  or  solid  modelling  of  the  milling  process,  in 
order  to  obtain  the  cutting  parameters.  Hence,  a reliable 
threshold  for  cutter  breakage  detection  can  be  established 
using  the  estimated  cutting  parameters  under  varying 
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cutting  conditions.  In  addition,  it  may  be  advantageous  to 
integrate  information  from  several  sensors  so  that  the  best 
features  of  all  the  sensors  are  utilized.  A multi-sensor 
scheme  using  artificial  intelligence  technologies  through  a 
knowledge-based  (analytical  and  experimental)  system  is 
also  considered.  Then,  the  tool  breakage  sensor  will  be 
fully  developed  to  achieve  the  final  goal  of  unmanned  milling. 


APPENDIX  A 


LISTING  OF  SUBROUTINES  FOR  MODELLING  OF  MILLING  PROCESS 


C SIMULATION  PROGRAM  FOR  MODELLING  OF  MILLING  PROCESS 
C 

C BY  YEONG-SHIN  TARNG 

C 

C 

C MILLIE. LIB  IS  A LIBRARY  FILE  CONTAINING  23  SUBROUTINES 

C 

C 


c 

c 

1) 

ANGLE 

— 

CALCULATE  ANGULAR  POSITIONS  OF  THE  TEETH  IN  THE 
INITIAL  CONDITION.  P.168. 

c 

c 

2) 

ROTAT 

— 

CALCULATE  INCREMENTAL  POSITIONS  OF  THE  TEETH  IN 
THE  CUTTING  OPERATION.  P.168. 

c 

c 

3) 

FORTR 

CALCULATE  TANGENTIAL  AND  RADIAL  FORCES  ON  EACH 
TOOTH.  (RIGID  FORCE  MODEL)  P.169. 

c 

c 

4) 

FXYS 

DECOMPOSE  TANGENTIAL  AND  RADIAL  FORCES  INTO  X 
AND  Y DIRECTIONS.  P.169. 

c 

c 

5) 

SUM 

SUM  THE  SIGNALS  IN  THE  X AND  Y DIRECTIONS. 
P. 170. 

c 

6) 

RESXY 

- 

CALCULATE  THE  RESULTANT  AVERAGE  SIGNAL.  P.170. 

c 

c 

7) 

BREAK 

— 

SIMULATE  THE  CHANGES  OF  TANGENTIAL  AND  RADIAL 
FORCES  DUE  TO  TOOL  BREAKAGE.  P.170. 

c 

c 

8) 

VI  BRA 

— 

SIMULATE  THE  VIBRATION  SIGNALS  IN  THE  X AND  Y 
DIRECTIONS.  P.171. 

c 

c 

9) 

THSEl 

CALCULATE  STARTING  AND  EXIT  ANGLES  OF  CUT  IN 
UP/DOWN  MILLING.  P.172. 

c 

c 

10) 

THSE2 

CALCULATE  STARTING  AND  EXIT  ANGLES  OF  CUT  IN 
UP/DOWN  MILLING.  P.173. 

c 

c 

11) 

THSE4 

* 

CALCULATE  STARTING  AND  EXIT  ANGLES  OF  CUT  IN 
UP/DOWN  MILLING.  P.174. 

c 

c 

12) 

THSE8 

CALCULATE  STARTING  AND  EXIT  ANGLES  OF  CUT  IN 
UP/DOWN  MILLING.  P.175. 

c 

c 

13) 

LOFIL 

* 

SIMULATE  THE  SIGNALS  PASSING  THROUGH  A 1ST  ORDER 
LOW-PASS  FILTER.  P.177. 

c 

c 

14) 

HIFIL 

— 

SIMULATE  THE  SIGNALS  PASSING  THROUGH  A 2ND  ORDER 
HIGH-PASS  FILTER.  P.177. 

c 

c 

15) 

FORTR 

CALCULATE  TANGENTIAL  AND  RADIAL  FORCES  ON  EACH 
TOOTH.  (END  MILLING  MODEL)  P.177. 

c 

c 

16) 

TRANS 

CALCULATE  INITIAL  ANGULAR  POSITIONS  OF  THE 
FLUTES.  P.178. 

c 

c 

17) 

HELIX 

— 

CALCULATE  INCREMENTAL  POSITIONS  OF  THE  FLUTES 
IN  END  MILLING.  P.178. 

c 

c 

18) 

UNROL 

— 

CALCULATE  THE  COEFFICIENTS  OF  THE  UNROLLED  END 
MILL  CUTTER.  P.179. 

c 

c 

19) 

FORTR 

— 

CALCULATE  TANGENTIAL  AND  RADIAL  FORCES  ON  EACH 
TOOTH.  (RUN-OUT  MODEL)  P.179. 

c 

20) 

THROW 

- 

INPUT  THE  ARRAY  OF  RADIAL  THROW.  P.180. 

c 

c 

21) 

FORTR 

— 

CALCULATE  TANGENTIAL  AND  RADIAL  FORCES  ON  EACH 
TOOTH.  (REGENERATIVE  FORCE  MODEL)  P.180. 

c 

c 

22) 

TEETH 

CALCULATE  THE  FEED  PER  TOOTH  FOR  EACH  TOOTH. 
(REGENERATIVE  FORCE  MODEL)  P.181. 

c 

c 

23) 

DEFLE 

SIMULATE  THE  VIBRATION  OF  CUTTER.  P.182. 
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C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


10 


C 


20 


C 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 


10 


c 


20 


SUBROUTINE  ANGLE ( THETA , THP , THO , NTEETH ) 

MILLING  SUBROUTINE 

ANGULAR  POSITIONS  OF  THE  TEETH 

THETA  - ANGULAR  POSITIONS  OF  THE  TEETH 
THP  - ANGLE  OF  THE  TOOTH  PERIOD 
THO  - INITIAL  SAMPLING  ANGLE 
NTEETH  - NUMBER  OF  THE  TEETH 

DIMENSION  THETA(l) 

DO  10  1=1, NTEETH 
THETA ( I ) = ( I - 1 ) *THP+THO 

DO  20  1=1, NTEETH 

IF  (THETA(I)  .GE.  360.)  THEN 
THETA ( I ) =THETA ( I ) -3  6 0 . 0 
ELSE 

IF  (THETA(I)  .LT.  0.)  THEN 
THETA(I)=THETA(I) +360.0 
ENDIF 
ENDIF 
CONTINUE 

RETURN 

END 


SUBROUTINE  ROTAT (THETA , THADD , THALL, NTEETH) 
MILLING  SUBROUTINE 

INCREMENTAL  POSITIONS  FOR  EACH  TOOTH 

THETA  - ANGULAR  POSITIONS  OF  THE  TEETH 
THADD  - INCREMENTAL  ANGLE  IN  THE  SIMULATION 
THALL  - TOTAL  ROTATION  ANGLE  IN  THE  CUTTING 
NTEETH  - NUMBER  OF  THE  TEETH 

DIMENSION  THETA (1) 

DO  10  1=1, NTEETH 
THETA ( I ) =THETA ( I ) +THADD 

DO  20  1=1, NTEETH 
IF  (THETA(I)  .GE.  360.)  THEN 
THETA ( I ) =THETA ( I ) -3  60 . 

ENDIF 

CONTINUE 


C 


o o 
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THALL=THALL+THADD 

RETURN 

END 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 


10 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


SUBROUTINE  FORTR ( FT , FR , THETA , THS , THE , NTEETH ) 

MILLING  SUBROUTINE  (RIGID  FORCE  MODEL) 

TANGENTIAL  & RADIAL  FORCES  ON  EACH  TOOTH 

FT  - ARRAY  OF  THE  TANGENTIAL  FORCES 

FR  - ARRAY  OF  THE  RADIAL  FORCES 

THETA  - ANGULAR  POSITIONS  OF  THE  TEETH 
THS  - STARTING  ANGLE  OF  CUT 
THE  - EXIT  ANGLE  OF  CUT 
NTEETH  - NUMBER  OF  THE  TEETH 

COMMON  /COEFF/  AXL, FEED, FK, RATIO, RAD 

DIMENSION  FT(1) ,FR(1) ,THETA(1) 

DO  10  1=1, NTEETH 

IF  (THETA(I)  .GE.  THS  .AND.  THETA(I)  .LE.  THE)  THEN 
FT ( I ) =FK*AXL*FEED*SIND (THETA ( I ) ) 

FR(I)=FT(I) *RATIO 
ELSE 

FT(I)=0.0 

FR(I)=0.0 

ENDIF 

CONTINUE 

RETURN 

END 


SUBROUTINE  FXYS ( FX , FY , FT , FR , THETA , NTEETH ) 
MILLING  SUBROUTINE 

DECOMPOSE  FORCES  INTO  THE  X & Y DIRECTIONS 

FX  - CUTTING  FORCE  IN  THE  X DIRECTION 

FY  - CUTTING  FORCE  IN  THE  Y DIRECTION 

FT  - ARRAY  OF  THE  TANGENTIAL  FORCES 

FR  - ARRAY  OF  THE  RADIAL  FORCES 

THETA  - ANGULAR  POSITIONS  OF  THE  TEETH 
NTEETH  - NUMBER  OF  THE  TEETH 


DIMENSION  FT(1) ,FR(1) ,THETA(1) 
C 

FSX=0 . 0 


ooo  oo  o o ooonoooono  on  o ooooooooon  no 
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FSY=0. 0 
C 

DO  10  I=1,NTEETH 

FSX=FSX+FR(I) *SIND(THETA(I) ) +FT ( I ) *C0SD (THETA ( I ) ) 
10  FSY=FSY-FR(I) *C0SD (THETA ( I ) )+FT(I) *SIND (THETA ( I ) ) 

C 

FX=FSX 

FY=FSY 

C 

RETURN 

END 


SUBROUTINE  SUM ( XCOM , YCOM , XAVG , YAVG) 

MILLING  SUBROUTINE 

SUM  THE  SIGNALS  IN  THE  X & Y DIRECTIONS 

XCOM  - COMPONENT  IN  THE  X DIRECTION 

YCOM  - COMPONENT  IN  THE  Y DIRECTION 

XAVG  - SUMMATION  OF  THE  X DIRECTION  COMPONENTS 

YAVG  - SUMMATION  OF  THE  Y DIRECTION  COMPONENTS 

XAVG=XAVG+DIMX 

YAVG=YAVG+DIMY 

RETURN 

END 


SUBROUTINE  RESXY (RES , XCOM, YCOM, NUM) 

MILLING  SUBROUTINE 

RESULTANT  AVERAGE  SIGNAL 

RES  - RESULTANT  SIGNAL 
XCOM  - SIGNAL  IN  THE  X DIRECTION 
YCOM  - SIGNAL  IN  THE  Y DIRECTION 
NUM  - NUMBER  OF  SAMPLES 


RNUM=NUM 

RES=(XCOM**2.+YCOM**2. ) **0.5/RNUM 

RETURN 

END 


SUBROUTINE  BREAK ( FT , FR , NTEETH , NBR , NP , PER , THETA) 
MILLING  SUBROUTINE 


u u 
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C CHIPPED  OR  TOTALLY  BROKEN  ON  THE  CUTTING  EDGE 

C 

C FT  - ARRAY  OF  THE  TANGENTIAL  FORCES 

C FR  - ARRAY  OF  THE  RADIAL  FORCES 

C NTEETH  - NUMBER  OF  THE  TEETH 

C NBR  - TOOTH  PERIOD  AS  THE  OCCURRENCE  OF  TOOL  BREAKAGE 

C NP  - CURRENT  TOOTH  PERIOD 

C PER  - PERCENTAGE  OF  THE  BROKEN  CUTTING  EDGE 

C THETA  - ANGULAR  POSITIONS  THE  TEETH 

C 

COMMON  /COEFF/  AXL, FEED, FK, RATIO, RAD 
C 

DIMENSION  FT(1) ,FR(1) ,THETA(1) 

C 

PERCENT=PER/100 . 

ERROR=FEED* PERCENT 
CHIPT=FEED*SIND (THETA (1) ) 

C 

IF  (NP  .GE.  NBR)  THEN 

IF  (CHIPT  .LE.  ERROR)  THEN 
FT(1)=0.0 
FR(1)=0.0 

FT (NTEETH) =2 . *FT (NTEETH) 

FR (NTEETH) =2 . *FR (NTEETH) 

ELSE 

FT (1)=(1. -PERCENT) *FT(1) 

FR(1)=(1. -PERCENT) *FR(1) 

FT ( NTEETH ) = (!.+ PERCENT ) * FT ( NTEETH ) 

FR ( NTEETH ) = (!.+ PERCENT ) * FR ( NTEETH ) 

ENDIF 

ENDIF 

C 

RETURN 

END 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


SUBROUTINE  VIBRA(FX,FY) 

MILLING  SUBROUTINE 

VIBRATIONS  IN  THE  X & Y DIRECTIONS 


FX 

FY 

DMASS 

ETA 

WN 

DDX,DDY 

DX,DY 

X,Y 

DT 


CUTTING  FORCE  IN  THE  X DIRECTION 
CUTTING  FORCE  IN  THE  Y DIRECTION 
MASS  PARAMETER 
DAMPING  RATIO 
NATURAL  FREQUENCY 

ACCELERATIONS  IN  X AND  Y DIRECTIONS 
VELOCITIES  IN  X AND  Y DIRECTIONS 
DISPLACEMENTS  IN  X AND  Y DIRECTIONS 
TIME  INTERVAL 


C 


COMMON  /VIBPAR/  WN, ETA, DMASS , DX, DY, DT, X, Y, TS 


n o 
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DDX=FX/DMASS-2 . *ETA*WN*DX- (WN**2 . ) *X 
DDY=FY/DMASS-2 . *ETA*WN*DY- (WN**2 . ) *Y 
C 

DX=DX+DDX*DT 

DY=DY+DDY*DT 

C 

X=X+DX*DT 

Y=Y+DY*DT 

C 

RETURN 

END 


SUBROUTINE  THSE 1 ( THS , THE , THP , THALL , WSL , DELY ) 

C 

C 1/1  IMMERSION  MILLING 

C 

C STARTING  & EXIT  ANGLES  FOR  CUTTING 

C 

C THS  - STARTING  ANGLE  OF  CUT 

C THE  - EXIT  ANGLE  OF  CUT 

C THP  - ANGLE  OF  THE  TOOTH  PERIOD 

C THALL  - TOTAL  ROTATION  ANGLE  IN  THE  CUTTING 

C WSL  - LENGTH  OF  THE  WORKPIECE 

C DELY  - DISTANCE  BETWEEN  THE  CUTTER  AND  WORKPIECE 

C 

COMMON  /COEFF/  AXL, FEED, FK, RATIO, RAD 
C 

DIMENS ION  THS ( 2 ) , THE ( 2 ) 

C 

FLEN=FEED* (THALL/THP) -DELY 
ELEN=WSL+RAD 
C 

DO  10  1=1,2 

THS (I) =0.0 
10  THE (I) =0.0 

C 

IF  (FLEN  .LT.  0.0)  THEN 
THS (1) =0.0 
THE(1)=0.0 
ELSE 

IF  (FLEN  .LT.  RAD)  THEN 

THCUT=ACOSD ( 1 . - (FLEN/ RAD) ) 

THS(1)=90.-THCUT 
THE ( 1 ) =9  0 . +THCUT 
ELSE 

IF  (FLEN  .LT.  WSL)  THEN 
THS (1) =0.0 
THE (1) =180.0 
ELSE 

THCUT=ACOSD ( 1 . - ( FLEN-WSL) /RAD) 

THS (1) =0.0 

THE(1)=90.-THCUT 

THS(2)=90.+THCUT 
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THE(2)=180.0 

ENDIF 

ENDIF 

ENDIF 

C 

RETURN 

END 


SUBROUTINE  THSE2 (THS , THE , THP, THALL, WSL, IUD,DELY) 

C 

C 1/2  IMMERSION  UP  AND  DOWN  MILLINGS 

C 

C STARTING  & EXIT  ANGLES  FOR  CUTTING 

C 

C THS  - STARTING  ANGLE  OF  CUT 

C THE  - EXIT  ANGLE  OF  CUT 

C THP  - ANGLE  OF  THE  TOOTH  PERIOD 

C THALL  - TOTAL  ROTATION  ANGLE  IN  THE  CUTTING 

C WSL  - LENGTH  OF  THE  WORKPIECE 

C lUD  - UP  OR  DOWN  MILLING 

C DELY  - DISTANCE  BETWEEN  THE  CUTTER  AND  WORKPIECE 

C 

COMMON  /COEFF/  AXL, FEED, FK, RATIO , RAD 
C 

FLEN=FEED* (THALL/THP) -DELY 
ELEN=WSL+RAD 
C 

IF  (FLEN  .LT.  0.0)  THEN 
THS=0 . 0 
THE=0 . 0 
ELSE 

IF  (lUD  .EQ.  1)  THEN 

IF  (FLEN  .LT.  RAD)  THEN 

THS=ASIND(1.-(FLEN/RAD) ) 

THE=90. 

ELSE 

IF  (FLEN  .LT.  WSL)  THEN 
THS=0 . 0 
THE=90.0 
ELSE 

IF  (FLEN  .LT.  ELEN)  THEN 
THS=0 . 0 

THE=90 . O-ACOSD ( 1 . “ ( FLEN-WSL) /RAD) 
ELSE 

THS=0 . 0 
THE=0 . 0 
ENDIF 
ENDIF 
ENDIF 
C 

ELSE 

C 

IF  (FLEN  .LT.  RAD)  THEN 
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C 

C 


THS=90. 0 

THE=180.0-ASIND(1.-(FLEN/RAD) ) 

ELSE 

IF  (FLEN  .LT.  WSL)  THEN 
THS=90.0 
THE=180.0 
ELSE 

IF  (FLEN  .LT.  ELEN)  THEN 

THS=90 . O+ACOSD ( 1 . - (FLEN-WSL) /RAD) 
THE=180. 0 
ELSE 

THS=0.0 
THE=0 . 0 
ENDIF 
ENDIF 
ENDIF 
ENDIF 

ENDIF 

RETURN 

END 


SUBROUTINE  THSE4 (THS , THE , THP, THALL, WSL, lUD, DELY) 

C 

C 1/4  IMMERSION  UP  AND  DOWN  MILLINGS 

C 

C STARTING  & EXIT  ANGLES  FOR  CUTTING 

C 

C THS  - STARTING  ANGLE  OF  CUT 

C THE  - EXIT  ANGLE  OF  CUT 

C THP  - ANGLE  OF  THE  TOOTH  PERIOD 

C THALL  - TOTAL  ROTATION  ANGLE  IN  THE  CUTTING 

C WSL  - LENGTH  OF  THE  WORKPIECE 

C lUD  - UP  OR  DOWN  MILLING 

C DELY  - DISTANCE  BETWEEN  THE  CUTTER  AND  WORKPIECE 

C 

COMMON  /COEFF/  AXL, FEED, FK, RATIO, RAD 
C 

FLEN=FEED* (THALL/THP) -DELY 
SLEN=(2.0-SQRT(3.0) )/2.0*RAD 
BLEN=WSL+SLEN 
ELEN=WSL+RAD 
C 

IF  (FLEN  .LT.  SLEN)  THEN 
THS=0 . 0 
THE=0 . 0 
ELSE 

IF  (lUD  .EQ.  1)  THEN 

IF  (FLEN  .LT.  RAD)  THEN 
THS=ASIND ( 1 . -FLEN/RAD) 

THE=ACOSD(l./2. ) 

ELSE 
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C 

C 


C 

C 


IF  (FLEN  .LT.  BLEN)  THEN 
THS=0.0 

THE=AC0SD(1./2. ) 

ELSE 

IF  (FLEN  .LT.  ELEN)  THEN 
THS=0 . 0 

THE=ASIND ( 1 . - (FLEN-WSL) /RAD) 

ELSE 

THS=0.0 

THE=0.0 

ENDIF 

ENDIF 

ENDIF 

ELSE 

IF  (FLEN  .LT.  RAD)  THEN 

THS=180. 0-ACOSD(l./2 . ) 

THE=180 . 0-ASIND ( 1 . -FLEN/RAD) 

ELSE 

IF  (FLEN  .LT.  BLEN)  THEN 

THS=180. 0-ACOSD(l./2 . ) 

THE=180.0 

ELSE 

IF  (FLEN  .LT.  ELEN)  THEN 

THS=180 . 0-ASIND ( 1 . “ ( FLEN-WSL) /RAD) 
THE=180.0 
ELSE 

THS=0 . 0 
THE=0 . 0 
ENDIF 
ENDIF 
ENDIF 
ENDIF 

ENDIF 

RETURN 

END 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


SUBROUTINE  THSE8 ( THS , THE , THP , THALL , WSL , lUD , DELY ) 

1/8  IMMERSION  UP  AND  DOWN  MILLINGS 

STARTING  & EXIT  ANGLES  FOR  CUTTING 

THS  - STARTING  ANGLE  OF  CUT 

THE  - EXIT  ANGLE  OF  CUT 

THP  - ANGLE  OF  THE  TOOTH  PERIOD 

THALL  — TOTAL  ROTATION  ANGLE  IN  THE  CUTTING 

WSL  - LENGTH  OF  THE  WORKPIECE 

lUD  - UP  OR  DOWN  MILLING 

DELY  - DISTANCE  BETWEEN  THE  CUTTER  AND  WORKPIECE 


COMMON  /COEFF/  AXL, FEED, FK, RATIO , RAD 

FLEN=FEED* (THALL/THP) -DELY 
SLEN=(4.-SQRT(7. ) )/4.*RAD 
BLEN=WSL+SLEN 
ELEN=WSL+RAD 

IF  (FLEN  .LT.  SEEN)  THEN 
THS=0.0 
THE=0 . 0 
ELSE 

IF  (lUD  .EQ.  1)  THEN 

IF  (FLEN  .LT.  RAD)  THEN 
THS=ASIND ( 1 . -FLEN/RAD) 

THE=ACOSD(3./4.) 

ELSE 

IF  (FLEN  .LT.  BLEN)  THEN 
THS=0 . 0 

THE=ACOSD(3./4. ) 

ELSE 

IF  (FLEN  .LT.  ELEN)  THEN 
THS=0.0 

THE=ASIND ( 1 . - ( FLEN-WSL) /RAD) 

ELSE 

THS=0 . 0 
THE=0 . 0 
ENDIF 
ENDIF 
ENDIF 

ELSE 

IF  (FLEN  .LT.  RAD)  THEN 
THS=180.0-ACOSD(3./4. ) 

THE=180 . 0-ASIND ( 1 . -FLEN/RAD) 

ELSE 

IF  (FLEN  .LT.  BLEN)  THEN 
THS=180.0-ACOSD(3./4. ) 

THE=180.0 

ELSE 

IF  (FLEN  .LT.  ELEN)  THEN 

THS=180 . 0-ASIND ( 1 . - ( FLEN-WSL) /RAD) 
THE=180. 0 
ELSE 

THS=0 . 0 
THE=0 . 0 
ENDIF 
ENDIF 
ENDIF 
ENDIF 


ENDIF 
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return 

END 
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SUBROUTINE  LOFIL (XI , YO , XI , Y1 , COEFl , COEF2 ) 
MILLING  SUBROUTINE 

1ST  ORDER  BUTTERWORTH  LOWPASS  FILTER 

XI  - INPUT  SIGNAL 
YO  - OUTPUT  SIGNAL 

YO=COEFl* Y1+COEF2  *XI 

X1=XI 

Yl=YO 

RETURN 

END 


SUBROUTINE  HIFIL (XI , YO , X2 , XI , Y2 , Yl) 

C 

C MILLING  SUBROUTINE 

C 

C 2ND  ORDER  BUTTERWORTH  HIPASS  FILTER 

C 

C XI  - INPUT  SIGNAL 

C YO  - OUTPUT  SIGNAL 

C 

COMMON  /HIPAS/  COEFl , COEF2 , COEF3 
C 

YO= (X2-2 . *X1+XI-C0EF1*Y2-C0EF2*Y1) /COEF3 
C 

X2=X1 

X1=XI 

C 

Y2=Y1 

Yl=YO 

C 

RETURN 

END 

C 

C 

SUBROUTINE  FORTR ( FT , FR , THETA , THS , THE , NTEETH , DAXL) 
C 

C MILLING  SUBROUTINE  (END  MILLING  MODEL) 

C 

C TANGENTIAL  & RADIAL  FORCES  ON  EACH  TOOTH 

C 

C FT  - ARRAY  OF  THE  TANGENTIAL  FORCES 

C FR  - ARRAY  OF  THE  RADIAL  FORCES 

C THETA  - ANGULAR  POSITIONS  THE  TEETH 
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C THS  - STARTING  ANGLE  OF  CUT 

C THE  - EXIT  ANGLE  OF  CUT 

C NTEETH  - NUMBER  OF  THE  TEETH 

C DAXL  - INCREMENTAL  AXIAL  DEPTH  OF  CUT 

C 

COMMON  /COEFF/  AXL, FEED, FK, RATIO , RAD 
C 

DIMENSION  FT(1) ,FR(1) ,THETA(1) 

C 

DO  10  1=1, NTEETH 

IF  (THETA(I)  .GE.  THS  .AND.  THETA(I)  . LE . THE)  THEN 
CHIPT=FEED*SIND(THETA(I) ) 

FT ( I ) =FK*DAXL*CHIPT 
FR(I)=FT(I) *RATIO 
ELSE 

FT (I) =0.0 
FR(I)=0.0 
ENDIF 
10  CONTINUE 
C 

RETURN 

END 

C 

C 

SUBROUTINE  TRANS ( PHI , THETA , NTEETH) 

C 

C MILLING  SUBROUTINE  (END  MILLING) 

C 

C ANGULAR  POSITIONS  OF  END  MILL  CUTTING  EDGES 

C 

C PHI  - ANGULAR  POSITIONS  OF  THE  CUTTING  EDGES 

C THETA  - STARTING  ANGULAR  POSITIONS  OF  THE  TEETH 

C NTEETH  - NUMBER  OF  THE  TEETH 

C 

DIMENSION  PHI ( 1 ) , THETA ( 1 ) 

C 

DO  10  1=1, NTEETH 
PHI(I)=THETA(I) 

10  CONTINUE 
C 

RETURN 

END 


SUBROUTINE  HELIX (PHI , NTEETH , DPS I) 

C 

C MILLING  SUBROUTINE  (END  MILLING) 

C INCREMENTAL  ANGULAR  POSITIONS  OF  END  MILL  CUTTING  EDGES 

C 

C PHI  - ANGULAR  POSITIONS  OF  THE  CUTTING  EDGES 

C NTEETH  - NUMBER  OF  THE  TEETH 

C DPS I - INCREMENTAL  ENGAGEMENT  ANGLE 

C 
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DIMENSION  PHI(l) 

C 

DO  10  I=1,NTEETH 

PHI(I)=PHI(I)-DPSI 
IF  (PHI(I)  .LT.  0.0)  THEN 
PHI(I)=PHI(I)+360. 
ENDIF 
10  CONTINUE 
C 

RETURN 

END 


SUBROUTINE  UNROL ( BETA, PSI , DPSI , DAXL, NSTEP) 

MILLING  SUBROUTINE  (END  MILLING) 

COEFFICIENTS  OF  END  MILL  CUTTER 

BETA  - HELIX  ANGLE 

PSI  - ENGAGEMENT  ANGLE 

DPSI  - INCREMENTAL  ENGAGEMENT  ANGLE 

DAXL  - INCREMENTAL  AXIAL  DEPTH  OF  CUT 

NSTEP  - NUMBER  OF  STEPS  IN  THE  SIMULATION 

COMMON  /COEFF/  AXL, FEED, FK, RATIO , RAD 

DIMENSION  PHI ( 1) , THETA (1) 

PSIRAD=AXL*TAND ( BETA) /RAD 
PSI=PSIRAD* (180. /3. 14159) 

DPSI=PSI/NSTEP 

DAXL=AXL/NSTEP 

RETURN 

END 


SUBROUTINE  FORTR ( FT , FR , THETA , THS , THE , NTEETH , ERROR) 
C 

C MILLING  SUBROUTINE  (RUN-OUT  MODEL) 

C 

C TANGENTIAL  & RADIAL  FORCES  ON  EACH  TOOTH 

C 

C FT  - ARRAY  OF  THE  TANGENTIAL  FORCES 

C FR  - ARRAY  OF  THE  RADIAL  FORCES 

C THETA  - ANGULAR  POSITIONS  OF  THE  TEETH 

C THS  - STARTING  ANGLE  OF  CUT 

C THE  - EXIT  ANGLE  OF  CUT 

C NTEETH  - NUMBER  OF  THE  TEETH 

C ERROR  - ARRAY  OF  THE  RADIAL  THROW 

C 

COMMON  /COEFF/  AXL, FEED, FK, RATIO , RAD 
C 
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DIMENSION  FT(1) ,FR(1) ,THETA(1) , ERROR ( -16 ; 16 ) 

C 

DO  30  I=1,NTEETH 

IF  (THETA(I)  .GE.  THS  .AND.  THETA(I)  . LE . THE)  THEN 
CHIP=FEED*SIND(THETA(I) ) 

CHIPT=CHIP+ERROR(I) -ERROR(I-l) 

RESUM=0 . 0 

IF  (CHIPT  .GT.  0.0)  THEN 
DO  10  J=1,NTEETH-1 

RES IDAL=CHI P+ERROR ( I - J ) -ERROR ( I - J- 1 ) 

IF  (RESIDAL  .LT.  0.0)  THEN 
RESUM=RESUM+RESIDAL 
ELSE 

GOTO  20 
ENDIF 

10  CONTINUE 

ELSE 

CHIPT=0. 0 
ENDIF 

20  FT(I)=FK*AXL* (CHIPT+RESUM) 

FR(I)=FT(I) *RATIO 
ELSE 

FT(I)=0.0 
FR(I)=0.0 
ENDIF 
3 0 CONTINUE 
C 

RETURN 

END 


SUBROUTINE  THROW ( ERROR , NTEETH) 

MILLING  SUBROUTINE 

RESET  THE  ARRAY  OF  RADIAL  THROW 

ERROR  - ARRAY  OF  THE  RADIAL  THROW 
NTEETH  - NUMBER  OF  THE  TEETH 

DIMENSION  ERROR(-16: 16) 

DO  10  1=1, NTEETH 
1 0 ERROR ( -NTEETH+ I ) =ERROR ( I ) 

RETURN 

END 


SUBROUTINE  FORTR ( FT , FR , THETA , THS , THE , NTEETH , FPT , Z BACK ) 
MILLING  SUBROUTINE  (REGENERATIVE  FORCE  MODEL) 
TANGENTIAL  & RADIAL  FORCES  ON  EACH  TOOTH 
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C 

C 

C 

C 

C 

C 

C 

c 

c 

c 


c 

c 

c 


c 


c 

10 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FT  - ARRAY  OF  THE  TANGENTIAL  FORCES 

FR  - ARRAY  OF  THE  RADIAL  FORCES 

THETA  - ANGULAR  POSITIONS  OF  THE  TEETH 

THS  - STARTING  ANGLE  OF  CUT 

THE  - EXIT  ANGLE  OF  CUT 

NTEETH  - NUMBER  OF  THE  TEETH 

FPT  - ARRAY  OF  THE  FEED  PER  TOOTH 

ZBACK  - SURFACE  PRODUCED  BY  THE  PRECEDING  CUT 

COMMON  /COEFF/  AXL, FEED, FK, RATIO , RAD 
COMMON  /SURWA/  Z , VIBX , VIBY , NZPOS , THADD 

DIMENSION  FT(1) ,FR(1) ,THETA(1) ,FPT(1) ,ZBACK(1) 

DO  10  1=1, NTEETH 

IF  (THETA(I)  .GE.  THS  .AND.  THETA(I)  .LE.  THE)  THEN 
NZPOS=THETA ( I ) /THADD 

Z =VIBX*SIND(THETA(I) )+VIBY*COSD(THETA(I) ) 

Z =Z*1000.0 

CHIPM=FPT(I) *SIND(THETA(I) ) 

CHIPT=CHIPM-ZBACK (NZPOS )+Z 
IF  (CHIPT  .LT.  0.0)  THEN 
CHIPT=0.0 

ZBACK (NZPOS )=ZBACK (NZPOS )-CHIPM 
ELSE 

ZBACK (NZPOS )=Z 
ENDIF 
ELSE 

CHIPT=0 . 0 
ENDIF 

FT ( I ) =FK*AXL*CHIPT 
FR(I)=FT(I) *RATIO 

CONTINUE 

RETURN 

END 


SUBROUTINE  TEETH ( FPT , NTEETH , NBR , NP , PER , THETA , ID) 

UNDAMAGED  OR  DAMAGED  CUTTER  (REGENERATIVE  FORCE  MODEL) 

FPT  - ARRAY  OF  THE  FEED  PER  TOOTH 
NTEETH  - NUMBER  OF  TEETH 

NBR  - TOOTH  PERIOD  AS  THE  OCCURRENCE  OF  TOOL  BREAKAGE 
NP  - CURRENT  TOOTH  PERIOD 

PER  - PERCENTAGE  OF  THE  BROKEN  CUTTING  EDGE 
THETA  - ANGULAR  POSITIONS  OF  THE  TEETH 
ID  - UNDAMAGED  OR  DAMAGED  CUTTER 


o o 


182 


COMMON  /COEFF/  AXL, FEED , FK, RATIO , RAD 
C 

DIMENSION  THETA ( 1) ,FPT(1) 

C 

DO  10  I=1,NTEETH 
10  FPT(I)=FEED 

C 

IF  (ID  .EQ.  2)  THEN 
PERCENT=PER/100 . 

ERROR=FEED* PERCENT 
CHIPT=FEED*SIND(THETA(1) ) 

IF  (NP  .GE.  NBR)  THEN 

IF  (CHIPT  .LE.  ERROR)  THEN 
FPT(l)  =0.0 

FPT (NTEETH) =2 . *FEED 
ELSE 

FPT(l)  = (1. -PERCENT) *FEED 

FPT (NTEETH) = ( 1 . +PERCENT) *FEED 
ENDIF 
ENDIF 
ENDIF 
C 

RETURN 

END 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 


C 


SUBROUTINE  DEFLE ( VIBX , VI BY , FX , FY , MODE , DT) 
MILLING  SUBROUTINE  (REGENERATIVE  FORCE  MODEL) 
VIBRATION  OF  THE  CUTTER 


VIBX  - VIBRATION  IN  THE  X DIRECTION 

VIBY  - VIBRATION  IN  THE  Y DIRECTION 

FX  - CUTTING  FORCE  IN  THE  X DIRECTION 

FY  - CUTTING  FORCE  IN  THE  Y DIRECTION 

MODE  - TOTAL  MODES  IN  THE  X AND  Y DIRECTIONS 

DT  - TIME  INTERVAL 

DMASS  - MASS  PARAMETER 

ETA  - DAMPING  RATIO 

WN  - NATURAL  FREQUENCY 

DDX,DDY  - ACCELERATIONS  IN  THE  X AND  Y DIRECTIONS 
DX,DY  - VELOCITIES  IN  THE  X AND  Y DIRECTIONS 
X,Y  - DISPLACEMENTS  IN  THE  X AND  Y DIRECTIONS 


COMMON  /SPINDL/  WN ( 10) , ETA( 10) , DMASS ( 10) 
COMMON  /EXCITX/  DDX ( 10) , DX ( 10) , X ( 10 ) 
COMMON  /EXCITY/  DDY ( 10) , DY ( 10) , Y ( 10) 


VIBX=0 . 0 
VIBY=0 . 0 


C 


DO  10  1=1, MODE 


DDX(I)=FX/DMASS(I) -2.*ETA(I) *WN(I) *DX(I) 
-(WN(I)**2.)*X(I) 

DDY(I) =FY/DMASS (I) -2 . *ETA(I) *WN(I) *DY(I) 
-(WN(I)**2.)*Y(I) 

DX(I)=DX(I)+DDX(I) *DT 
DY(I)=DY(I)+DDY(I) *DT 

X(I)=X(I)+DX(I) *DT 
Y(I)=Y(I)+DY(I) *DT 

VIBX=VIBX+X(I) 

VIBY=VIBY+Y(I) 


RETURN 


APPENDIX  B 


EXAMPLE  PROGRAM  FOR  SIMULATION  OF  DISPLACEMENT 
SIGNAL  OF  SPINDLE 


ooo  on  n o o ooooooo 


C 

c 

C 

C AN  EXAMPLE  SHOWS  HOW  TO  CALL  < MILLIE. LIB  > FOR 

C 

C SIMULATION  OF  DISPLACEMENT  SIGNAL  OF  SPINDLE 

C 
C 
C 

C CUTTING  PARAMETERS 

C 

COMMON  /COEFF/  AXL, FEED , FK, RATIO , RAD 
VIBRATION  PARAMETERS 

COMMON  /VIBPAR/  WN , ETA, DMASS , DX, DY , DT , X , Y , TS 

MAXIMUM  TEETH  ON  THE  CUTTER:  16 

DIMENSION  FT ( 16 ) , FR ( 16 ) , THETA ( 16 ) 

PRINT* NUMBER  OF  THE  TEETH  ON  THE  CUTTER' 
READ(*,*)  NTEETH 

PRINT*, '(1)  UP  OR  (2)  DOWN  MILLING  ->  (1,2)' 
READ(*,*)  lUD 

PRINT* ,' SPINDLE  SPEED(R.  P.  M.)' 

READ(*,*)  RPM 

THP  =3 60. /NTEETH 
DPS  =RPM*360./60. 

SPT  =THP/DPS 

PRINT*, 'SEC  PER  TOOTH  — > ',SPT 

PRINT*, 'SIMULATION  TIME  INTERVAL  (SEC)' 
READ(*,*)  DT 

PRINT* ,' SAMPLING  TIME  INTERVAL  (SEC)' 

READ(*,*)  TS 

PRINT* INITIAL  SAMPLING  ANGLE  (DEG)' 

READ(*,*)  THO 

CALCULATE  SIMULATION  TIME  COEFFICIENTS 

THADD=DT*DPS 
NPTS  =THP/THADD+0.5 
INTVL=TS/DT+0.5 
ICOUNT=INTVL-l 
NS AMP=NPTS/ INTVL 
C 

PRINT*, 'TOOL  BROKEN  ? ->  (Y/N) ' 
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C 


10 


C 


C 

C 

C 


C 

C 

C 


C 


READ (*,10)  TYPE 
FORMAT (Al) 

IF  (TYPE  .EQ.  'Y')  THEN 

PRINT* INPUT  FROM  WHICH  TOOTH  PERIOD 

READ ( * , * ) NBR 

PRINT*  ,*  PERCENTAGE  OF  MISSING  CUTTING  EDGE  (-s)' 
READ ( * , * ) PER 
ENDIF 

PRINT*, 'CLEARANCE  BETWEEN  THE  WORKPIECE  AND  CUTTER' 
READf*.*)  DELY 

PRINT* ,' TOTAL  TOOTH  PERIODS  IN  THE  SIMULATION 
READ(*,*)  NPMAX 

CUTTING  PARAMETERS 

AXL  =2.54 
FEED  =0.25 
FK  =1500.0 
RATIO=0. 5263 
RAD  =50.8 
WSL  =100.0 

VIBRATION  COEFFICIENTS 

DMASS=1.0 
ETA  =0.04 
FN  =600.0 
WN  =2 . *3 . 14159*FN 

TIME  =0.0 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


OPEN  A DATA  FILE  FOR  OUTPUT 
OPEN (UNIT=5 , FILE= ' STEP . DAT ' , STATUS= ' NEW ' ) 
CALCULATE  ANGULAR  POSITIONS  OF  THE  TEETH 
CALL  ANGLE ( THETA , THP , THO , NTEETH ) 

* stMUIATTON  of  the  CUTTING  PROCESS  * 

DO  20  1=1, NPMAX 
DO  20  J=1,NPTS 


CALCULATE  STARTING  AND  EXIT  ANGLES  OF  CUT 
IN  THE  1/2  IMMERSION  MILLING 

CALL  THSE2 (THS , THE , THP , THALL, WSL, lUD, DELY) 


C 
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C 

C 

C 

C 

C 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

20 

c 


CALCULATE  TANGENTIAL  AND  RADIAL  FORCES 
(FACE  MILLING,  RIGID  FORCE  MODEL) 

CALL  FORTR ( FT , FR , THETA , THS , THE , NTEETH ) 

MONITOR  TOOL  BREAKAGE 

IF  (TYPE  .EQ.  'Y')  THEN 

CALL  BREAK ( FT , FR , NTEETH , NBR , I , PER , THETA ) 

ENDIF 

DECOMPOSE  FORCES  INTO  X AND  Y DIRECTIONS 
CALL  FXYS ( FX , FY , FT , FR , THETA , NTEETH ) 

SIMULATE  VIBRATION  SIGNALS  IN  THE  X AND  Y DIRECTIONS 
CALL  VIBRA(FX,FY) 

ICOUNT=ICOUNT+l 

STORE  SIMULATION  DATA  INTO  THE  DATA  FILE 

IF  (ICOUNT  .EQ.  INTVL)  THEN 
ICOUNT=0 

WRITE(5,*)  TIME,X,Y 
ENDIF 

CALCULATE  INCREMENTAL  POSITIONS  OF  THE  TEETH 

CALL  ROTAT ( THETA , THADD , THALL , NTEETH ) 

TIME=TIME+DT 

STOP 

END 
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